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ABSTRACT 


A  review  of  the  literature  on  spall  fracture  and  dynamic  response  of 
materials  is  presented.  Current  theories  underlying  large  amplitude  wave  prop¬ 
agation  and  spall  fracture  under  conditions  of  uniaxial  strain  are  described  in 
detail.  Relevant  features  of  elastic,  plastic,  and  shock  wave  phenomena  ere 
reviewed,  and  the  effects  of  subsidiary  properties  on  the  material  response  are 
cons'.dered,  including  strain  rate,  temperature,  finite  strain,  work  hardening, 
and  the  Bauschinger  effect. 
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SPALL  FRACTURE  AND  DYNAMIC  RESPONSE  OF  MATERIALS 


by 

J,  H.  Oscarson 
X.  F.  Graff 

INTRODUCTION 

This  report  presents  a  review  of  the  literature  on  spall  fracture  and 
dynamic  response  of  materials.  Since  the  purpose  of  the  review  is  to  provide 
background  and  direction  for  an  experimental  and  analytical  spall  study  presently 
being  conducted  at  Battelle-Columbus  Laboratories  (BCL) ,  the  greatest  emphasis 
is  given  to  these  topics  most  closely  related  to  this  study.  With  this  con¬ 
straint,  most  of  the  subject  matter  will  be  devoted  to  a  detailed  description  of 
the  current  theories  underlying  large  amplitude  wave  propagation  and  spall 
fracture  under  conditions  of  uniaxial  strain. 

Fundamental  to  the  study  of  spallation  phenomena  is  determination 
of  stress  history  at  the  spall  plane.  However,  the  stress  cono4.>*.or.  in  the 
interior  of  the  material  cannot  be  measured  directly,  but  must  be  inferred  from 
measurements  of  free  surface  motions  of  the  target  material  and  solution  of 
the  stress  wave  propagation  problem.  It  is  necessary,  as  a  prelude  to  solution 
of  the  stress  wave  problem,  to  know  the  dynamic  properties  of  the  material — 
that  is,  the  constitutive  relation  of  the  material. 

The  purpose  of  the  appendices  is  to  review  the  relevant  features  of 
elastic,  plastic,  and  shock  wave  phenomena.  Prime  attention  is  directed  to  a 
review  of  elastic-plastic  stress  wave  theory  with  the  intent  of  relating  experi¬ 
mental  measurements  of  rear  surface  motion  to  the  gross  features  of  the  interior 
stress  history.  It  will  be  found,  however,  that  the  detailed  analysis  of  the 
stress  history  is  too  complex  for  simple  analysis  and  must  await  computer 
investigation. 

'Ihe  effects  of  a  number  of  subsidiary  properties  on  the  material 
response  will  also  be  considered.  These  will  include  the  effects  of  strain 
rate,  temperature,  finite  strain,  work  hardening,  arid  the  Bauschinger  effect. 

Some  assessment  of  these  effects  is  necessary  if  the  spall  study  data  are  to  be 
interpreted  unobscured  by  unknown  complexities  of  material  response. 


Unfortunately,  detailed  knowledge  of  dynamic  material  properties  is  scant,  so 
that  it  is  only  pcssible  to  estimate  the  consequences  of  a  number  of  the  material 
effects.  It  sboul.  be  noted  that  dynamic  material  properties  is  currently  the 
subject  of  intense  research. 

Ihe  early  time  response  of  a  plate  subjected  to  lateral  impulsive 
loading  takes  the  form  of  a  compressive  dilatational  stress  wave  propagating 
through  the  thickness.  If  the  induced  stresses  are  greater  than  the  elastic 
limit  of  the  material,  the  peak  stress  of  the  wave  will  be  attained  and  propa¬ 
gate  behind  a  higher  velocity  elastic  precursor  wave.  Ihus,  at  the  strain  rates 
of  interest  (10^  -  10^/sec),  the  unstressed  material  will  first  experience  an 
almost  instantaneous  rise  to  the  stress  corresponding  to  the  elastic  limit.  It 
will  remain  at  this  stress  for  some  finite  time,  and  then  experience  another 
rapid  rise  to  the  peak  stress.  Details  of  the  unloading  will  depend  in  large 
part  on  how  the  impulsive  load  was  generated,  i.e.,  contact  explosive  or  flying 
plate  impact. 

When  the  front  of  the  waveform  reaches  the  stress-free  rear  surface  of 
the  plate,  it  reflects  as  a  tensile  wave  and  propagates  into  the  compressed 
material  associated  with  that  part  of  the  wave  which  has  yet  to  be  reflected. 

If  the  total  duration  oi  the  incoming  compressive  pulse  is  sufficiently  short, 
there  will  be  a  plane  at  some  finite  distance  from,  and  parallel  to,  the  rear 
surface,  for  which  the  net  stress  first  becomes  tensile.  (Here,  net  stress 
refers  to  the  sum  of  the  compressive  st-ess  associated  with  the  incoming  wave 
and  the  tensile  stress  associated  with  the  reflected  wav?..)  If  the  stress 
and  pulse  duration  are  sufficiently  large,  the  material  may  fracture  (spall)  at 
or  near  this  plane.  Appendices  A,  B,  C,  D,  and  E  describe  all  the  events 
leading  up  to  conditions  sufficient  for  spall. 

"General  Physical  Features  of  Plastic  Flow  and  Fracture"  provides  a 
brief  and  very  elementary  review  of  the  physics  and  mechanics  of  flow  and  fracture 
from  a  primarily  micromechanical  point  of  view.  It  is  included  in  this  report  so 
that  those  readers  unfamiliar  with  these  concepts  might  acquire  some  background 
for  the  topics  discussed  in  subsequent  sections. 
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"Spall'*  is  a  detailed  statement  of  the  current  knowledge  regarding 
spall  fracture.  This  includes  a  discussion  of  the  factors  affecting  spall,  as 
well  as  the  means  by  which  these  factors  and  their  effects  are  used  to  predict 
spall. 

Hie  final  section  is  a  brief  summary  of  conclusions  drawn  from  the 
review  of  literature  on  spall. 

Appendices  A,  B,  C,  D,  and  E  present  a  detailed  discussion  of  large 
amplitude  wave  propagation  and  constitutive  relations  of  materials,  the  major 
emphasis  being  given  to  metals. 


DISCUSSION 


General  Physical  Features  of  Plastic 
Plow  and  Fracture 


Observations  generally  descriptive  of  material  behavior  in  standard 
tensile  or  compressive  tests  (i.e.,  quasi-static  loading  and  uniaxial  stress 
states)  are  sufficiently  numerous  and  well-developed  to  serve  as  a  background 
for  understanding  the  physical  nature  of  stress  wave  propagation  and  spallation 
in  uniaxial  strain.  Furthermore,  the  existence  and  operation  of  many  of  the 
basic  mechanisms  are  only  quantitefci vtly  modified  by  changes  in  loading  condi¬ 
tions.  Thus,  before  proceeding  to  detailed  discussions  on  material  response  and 
spall,  a  brief  review  of  the  essential  physical  features  of  flow  and  fracture 
appears  to  be  in  order. 

Riysics  and  Mechanics  of  Plastic  Flow 

Composition  of  Metals.  Metallic  elements  consist  of  atoms  in  an 
ordered  arrangement,  called  a  lattice,  at  equilibrium  under  their  mutual  forces. 
Typical  crystallographic  arrangements  ere  body-centered  cubic  (BCC)  and  face- 
centered  cubic  (FCC)  as  shown  in  Figure  1. 


FIGURE  1.  BODY-CENTERED  AND  FACE-CENTERED  LATTICE 
ARRANGEMENTS 
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These  figures  represent  only  a  small  portion  of  the  total  possible  lattice 
arrangements.  However,  over  70  percent  of  metals  and  alloys  crystallize  in 
the  above  or  in  a  hexagonal  system. 

A  pure  and  perfect  metallic  specimen  would  consist  of  an  uninterrupted 
repetition  of  the  particular  lattice  arrangement  of  the  metal.  However, 
typically,  as  a  liquid  metal  solidifies  from  a  melt,  a  variety  of  defects  are 
formed.  The  simplest  defect,  called  a  vaccncy,  occurs  when  an  atom  is  missing 
trom  the  lattice.  An  extra  atom  wedged  in  the  lattice  between  the  normal 
atomic  sites  may  also  occur;  this  is  termed  an  interstitial  atom.  Impurity 
atoms  may  replace  a  regular  lattice  atom  or  be  accommodated  within  the  lattice 
as  an  interstitial  impurity.  Such  defects,  involving  only  one  atom,  are  called 
point  defects.  Line  defects  are  the  next  level  of  lattice  imperfection.  The 
incomplete  translation  of  one  part  of  the  lattice  relative  to  the  other, 
designated  as  a  dislocation,  is  an  example. 

Surface  defects  within  the  crystal  may  also  exist.  Thus,  across  a 
surface  defect,  the  crystal  structure  may  undergo  a  change  of  orientation.  An 
example  would  be  the  boundary  which  occurs  when  the  crystallographic  structure 
forms  a  mirror  image  of  itself  as  shown  in  Figure  2. 

As  a  metal  solidifies  from  a  melt,  crystals  grow  from  many  nuclei  on 
cooling  so  that  the  resultant  solid  consists  of  many  grains  of  metal,  each  grain 
composed  of  the  crystal  lattice.  Since  the  grains  will  be  misoriented  relative 
to  one  anocher,  grain  boundaries  are  formed  where  some  compromise  of  lattice 


FIGURE  2.  TWIN  BOUNDARY  WITHIN  A  CRYSTAL 


r 
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position  occurs.  In  fact,  if  the  relative  misorientation  between  grains  is 
large,  the  boundary  may  be  considered  a  thin  transition  layer  of  high  disorder 
with  amorphous  rather  than  crystalline  properties.  The  situation  is  schemati¬ 
cally  illustrated  in  Figure  3. 


Microscopic  Deformation  of  Metals.  In  the  case  of  small  loads  applied 
to  a  perfect  crystal,  small  distortions  will  occur  in  the  crystal  lattice  as  the 
atoms  seek  new  equilibrium  positions  in  response  to  the  applied  forces.  If  the 
lattice  is  perfect,  enormous  forces  are  required  to  rupture  the  atomic  bonds — 
as  shown  by  tests  on  perfect  crystals.  Once  the  loads  are  removed,  the  atoms 
return  to  their  original  positions  and  the  material  thus  exhibits  "perfect 
elasticity".  The  case  for  polycrystals  is  similar  except  the  loads  necessary 
to  produce  some  irreversible  changes  in  the  lattice  structure  are  much  smaller; 
therefore,  polycrystals  are  also  capable  of  exhibiting  perfect  elasticity,  but 
to  a  lesser  extent. 

The  case  of  permanent  changes  in  the  lattice  structure  caused  by  loads 
greater  than  considered  above  will  now  be  considered.  The  fundamental  mechanisms 
of  permanent  or  plastic  deformation  within  a  crystal  are  the  generation  and  motion 
of  dislocations.  The  two  basic  dislocation  types  are  screw  dislocation  and 
edge  dislocation.  In  a  screw  dislocation,  relative  shifting  of  the  crystallo¬ 
graphic  planes  has  occurred,  but  this  has  not  propagated  through  the  entire 
crystal  and,  in  fact,  terminates  along  the  screw  dislocation  line.  The  edge 
dislocation,  which  is  simpler  to  visualize,  is  the  result  of  an  extra  half  plane 
of  atoms  being  "pushed"  into  the  lattice. 


FIGURE  3.  GRAINS  AND  GRAIN  BOUNDARIES  IN  A  POLYCRYSTAL 


In  either  type  of  dislocation,  a  plane  or  a  portion  of  a  plane  of 
atoms  has  glided  with  respect  to  an  adjacent  plane,  and  after  the  gliding 
motion,  the  geometric  pattern  of  the  crystal  has  again  come  back  into  registry. 

In  all  types  of  deformation  considered,  the  lattice  undergoes  some 
distortion  in  the  vicinity  of  the  crystal  fault;  consequently,  the  lattice 
possesses  additional  strain  energy  and  is  therefore  in  a  higher  energy  state 
than  without  the  fault.  Nevertheless,  the  essential  lattice  structure  is  only 
locally  distorted  by  the  presence  of  a  dislocation  and  the  identifying  features 
of  the  crystal  are,  in  the  main,  retained. 

Macroscopic  Deformation  of  Metals.  The  microscopic  deformation  of 
metals  manifests  itself  macroscopically  as  plastic  flow  and  deformation.  That 
is,  microscopic  dislocations  translate  as  plastic  deformations  of  the  crystal 
or  polycrystal.  Although  an  annealed  crystal  will  contain  a  large  number  of 
dislocations  per  unit  volume  (10^  to  10^  per  cubic  centimeter),  the  magnitudes 
of  macroscopically  observed  plastic  strains  are  not  solely  attributable  to  the 
gliding  of  these  existing  dislocations.  It  is  known  that  the  dislocation  density 
increases  with  the  start  of  plastic  deformation.  Although  the  description  of 
these  multiplication  mechanisms  is  beyond  the  scope  of  this  review,  both  the 
theoretical  basis  for  their  operation  and  X-ray  evidence  of  their  existence 
have  been  established.  Thus,  flow  takes  place  due  to  relative  motion  of  various 
dislocation  glide  planes.  The  number  of  operative  glide  planes  tends  to 
greatly  increase  due  to  multiplication  mechanisms  as  flow  progresses,  resulting 
in  macroscopic  plastic  deformation. 

Physics  and  Mechanics  of  Fracture 

Fracture  in  most  materials  requires  the  contribution  of  two  usually 
distinct  mechanisms;  crack  or  void  nucleation  and  crack  growth  or  propagation. 
Since  the  operation  and  characteristics  of  both  mechanisms  are  highly  dependent 
on  the  inelastic  behavior  modes  of  the  material,  the  discussion  to  follow  will 
consider  separately  the  fracture  of  perfectly-brittle  materials,  semibrittle 
materials,  and  ductile  materials.* 


*  Principal  sources  for  this  discussion  are  References  1,  2,  3,  and  4,  page  35. 


Perfectly-Brlttle  Materials 


Under  certain  conditions,  e.g.,  low  temperatures,  both  amorphous 
materials  and  some  metals  will  behave  in  a  semibrittle  or  ideally  brittle  fashion. 
(Fracture  of  materials  falling  in  the  semibrittle  category  are  discussed  in  a 
subsequent  section.)  Those  materials  belonging  to  the  latter  category,  which  have 
been  studied  most  extensively,  are  the  inorganic  glasses.  Materials  such  as 
ceramics  have  received  less  attention. 

Crack  Nucleation.  The  discrepancy  between  the  theoretical  cohesive 
strength  (5-10%  of  E*)  and  the  experimentally  determined  tensile  strength  is 
attributed  to  stress  concentrations  at  the  tips  of  randomly  distributed  microcracks 
on  the  surface  or  in  the  volume  of  the  material.  For  uniformly  stressed  material, 
the  most  severe  (longest)  crack  propagates  through  the  material  in  a  direction 
perpendicular  to  the  maximum  principal  normal  stress,  and  thus  causes  failure  at 
relatively  low  applied  stress.  While  these  microcracks  have  never  been  micro¬ 
scopically  observed,  order-of-magnitude  increases  in  the  strength  of  glass  specimens 
with  a  surface  layer  etched  off  have  been  obtained.  In  the  case  of  glass,  the 
severe  surface  cracks  which  cause  failure  are  attributed  to  mechanical  damage 
introduced  in  processing  or  handling. 

Crack  Propagation.  In  uniaxial  tension,  the  stress,  oq,  necessary  to 
propagate  an  elliptical  crack  of  length  2C  (with  major  axis  normal  to  ctq)  was 
derived  by  Griffith:  ^ 


where  a  «  specific  surface  energy.  When  this  equation  is  satisfied,  the  crack 
begins  to  propagate  across  the  plane  of  maximum  tension.  As  the  crack  elongates, 
the  stress  concentration  at  its  tip  increases,  thereby  accelerating  the  motion 


*Young's  Modulus. 
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of  the  tip  to  maximum  velocities  as  high  as  the  Rayleigh  wave  velocity.  As  long 
as  is  maintained  on  the  specimen,  the  crack  will  continue  to  move  at  steady- 
state  velocity  unless  it  hits  an  obstacle  (hard  region)  or  forks. 

Materials  with  crack  nucleation  and  propagation  characteristics  like 
those  described  above  exhibit  pronounced  size  effects  in  their  tensile  strengths, 
i.e.,  large  specimens  are  weaker  than  small.  This  is  due  to  the  fact  that  the 
probability  of  a  specimen  containing  a  severe  crack  increases  with  size. 

The  simple  considerations  introduced  above  have  been  extended  to  more 
elaborate  analyses;  these  include  biaxial  and  triaxial  stress  fields  and 
statistical  approaches  relating  flaw  densities  and  distributions  to  strength. 

Semlbrittle  Materials 


Some  body-centered  cubic  metals  (BCC) ,  such  as  low-carbon  steel  and 
tungsten,  and  glassy  materials,  such  as  inorganic  glasses  and  polymers,  will 
deform  plastically  or  viscously  at  high  temperatures  or  low  strain  rates,  but  at 
low  temperatures  or  high  strain  rates  they  will  fracture  in  a  semlbrittle  manner. 

In  both  crack  nucleation  and  propagation,  striking  differences  can  be  found 
between  semlbrittle  fracture  and  the  perfectly-brittle  fracture.  Another  factor 
which  contributes  to  the  ductile/ semlbrittle  fracture  mode  transition  of  these 
materials  is  the  presence  of  a  high  triaxial  tensile  stress.  In  the  poly¬ 
crystalline  materials,  the  fractures  are  usually  transcrystalline  (cleavage  of 
individual  crystals),  although  some  semlbrittle  fracture  is  of  the  intercrystalline 
variety  (brittle  separation  along  grain  boundaries). 

Crack  Nucleation.  The  main  difference  between  perfectly-brittle  and 
semlbrittle  fracture  in  the  glassy  materials  is  that,  in  the  latter  cnse,  some 
plastic  deformation  may  precede  or  accompany  failure.  The  mechanism  for  crack 
nucleation  is  the  same  as  that  discussed  earlier.  In  metals,  however,  cracks  are 
preceded  and  can  be  nucleated  by  inhomogeneous  plastic,  deformation  on  a 
microscopic  scale,  i.e.,  slip  and  twinning. 
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Crack  nucleation  by  slip  can  be  caused  by  an  edge  dislocation  pile-up 
at  a  strong  obstacle  such  as  a  high  or  medium  angle  grain  boundary  (see  Figure  4) 
or  a  pile-up  in  the  form  of  a  low  angle  boundary  (see  Figure  5).  In  both  cases, 
the  concentrated  stress  becomes  equal  to  the  theoretical  strength  of  the  material, 
and  as  long  as  the  stress  does  not  relax  too  rapidly,  the  cracks  will  continue  to 
grow.  An  interesting  aspect  of  the  grain  boundary  pile-up  mechanism  is  that  the 
crack  nucleation  stress  is  grain-size  dependent.  This  is  due  to  the  fact  that 
more  dislocations,  and  hence  greater  stress  concentrations,  will  occur  in  larger 
grained  materials.  Pile-ups  in  the  form  of  low-angle  boundaries  are  common  to 
hexagonal  metals  and  single  crystals  of  ionic  compounds,  such  as  magnesium  oxide. 
(See  Figure  5.) 

A  sketch  of  the  mechanism  of  crack  nucleation  by  the  interaction  of 
deformation  twins  is  shown  in  Figure  6.  This  mechanism,  found  in  BCC  crystals, 
requires  higher  stresses  than  those  needed  for  dislocation  pile-ups.  Only  a  few 
intersections  of  this  type  will  result  in  crack  nucleation,  since  the  stresses  are 
easily  relieved. 


Crack  Propagation.  Once  defects  (cracks)  have  been  introduced  by  the 

flow  of  crystalline  material,  the  occurrence  of  semibrittle  or  ductile  fracture 

will  be  governed  by  whether  or  not  an  elastic-plastic  crack  propagation  condition 

(1) 

is  satisfied.  Such  a  condition  was  given  by  Orowan  as 


2„E 


(2) 


o  ttC 

where  p  *  the  work  of  plastic  deformation  per  unit  increase  in  crack  area.  This 
is  Just  the  Griffith  equation  with  p  substituted  for  a.  In  general,  however,  p 
will  not  be  constant. 


Observed  from  a  microscopic  point  ot  view,  it  is  probable  that  a  flow- 
nucleated  crack  in  a  single  crystal  would  be  of  sufficient  length  to  satisfy  the 
Griffith  brittle  crack  propagation  criterion.  In  a  polycrystalline  material 
however,  a  cleavage  crack  will  usually  be  stopped  by  a  grain  boundary  between 
grains  with  greatly  differing  cleavage  planes.  Once  the  crack  is  stopped,  plastic 
deformation  at  its  tip  serves  to  blunt  it,  thus  creating  the  necessity  for 
additional  crack  propagation  energy.  For  this  reason,  Che  term  "p"  in  Orowan's 
equation  ' _  much  higher  than  the  corresponding  "or"  term  in  Griffith's  equation. 


FIGURE  4.  FORMATION  OF  A  CRACK  BY  FUSING  TOGETHER  OF 
DISLOCATIONS  AT  THE  TIP  OF  A  PILE-UP 


^  (b)  (OF  197) 

FIGURE  5.  FORMATION  OF  A  MICROCRACK  BY  SHEAR 
OF  PART  OF  A  TILT  BOUNDARY 


FIGURE  6.  CRACK  NUCLEATION  BY  THE  INTERSECTION 
OF  WO  DEFORMATION  WINS 
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As  was  the  case  for  crack  nucleatlon  in  polycrystalline  materials , 
grain  size  plays  an  important  role  in  crack  propagation.  The  improved  resistance 
to  semibrittle  fracture  of  fine-grained  materials  is  due  not  only  to  the  fact  that 
cracks  are  more  difficult  to  nucleate,  but  also  that  they  are  more  difficult  to 
propagate . 

Intercrystalline  semibrittle  fracture  is  nucleated  and  propagated  in 
the  grain  boundaries.  It  is  due  either  to  some  form  of  grain  boundary  ambritt..e- 
rnent  or  a  thin  film  of  a  brittle  phase  separating  at  the  grain  boundary. 

Ductile  Materials 


Materials  such  as  some  FFC  (face-centered  cubic)  metals  do  not  undergo 
a  ductile/semibrittle  fracture  mode  transition.  In  a  standard  tension  test,  these 
materials  develop  voids  or  ductile  cracks  under  the  influence  of  triaxial  tension 
in  the  center  of  the  necked  region.  These  voids  subsequently  grow  in  size  and 
number  and  coalesce  as  the  load  is  maintained  or  increased.  The  fine  structure 
of  the  fracture  surface  at  the  center  of  a  completely  separated  specimen  shows 
that  the  voids  grow  to  such  an  extent  that  the  material  between  them  ruptures 
(shear  failure)  rather  than  cleaves. 

Crack  Nucleation.  The  origin  of  the  voids  which  grow  to  cause  ductile 
fracture  is  attributed  to  either  stress  concentrations  in  the  grain  boundaries 
due  to  impurities  or  inclusions,  fracture  of  perfectly* brittle  elements  (such  as 
pearlite  in  steel),  or  excessive  localized  strain  concentrations  in  bands  of 
heavy  deformation. 

Crack  Propagation.  Crack  propagation,  or  more  descriptively,  void 
growth,  proceeds  according  to  the  manner  in  which  the  crack  was  nucleated.  Thus, 
if  the  cracks  were  initiated  by  failure  of  brittle  elements,  they  would  become 
blunted  and  turn  into  holes.  The  growth  and  coalescence  of  the  holes  would  depend 
on  the  amount  of  further  plastic  strain  and  the  magnitude  of  the  transverse  stress. 
On  the  other  hand,  if  the  cracks  were  formed  along  grain  boundaries  or  bands  of 
heavy  deformation,  the  cavity  acts  as  a  notch  to  concentrate  stress  in  planes 
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oriented  at  30  to  40  degrees  to  the  tensile  axis.  Large  amounts  of  plastic  flow 
occur  in  this  region,  forming  a  heavy  concentration  of  voids  in  the  shear  band, 
and  the  weakened  region  then  splits  open. 

Spall 

This  section  is  devoted  to  a  detailed  statement  of  the  current  know¬ 
ledge  of  spall  fracture.  For  reasons  of  readability  and  conceptual  clarity,  the 
results  of  the  various  Investigators  are  presented  in  an  integrated  (5)  rather 
than  sequential  (*>)  manner.  in  this  format,  topics  such  as  spall  models  and 
degrees  of  snail  are  discussed  separately.  Particular  emphasis  is  placed  on 
those  results  most  closely  related  to  spallation  in  uniaxial  strain  caused  by 
plate-slap  impact,  i.e,,  rectangular  pulses.  Furthermore,  the  closely  related 
and  at  least  equally  important  problem  of  determining  material  response  prior  to 
spall  will  receive  only  cursory  attention,  since  this  has  been  considered  in 
detail  in  Appendices  A  to  E. 

The  Degrees  of  Spall 

The  most  general  definition  of  spall  is  given  in  terms  of  its  cause; 
namely,  material  failure  due  to  the  interaction  of  two  or  more  rarefaction 
waves.  The  failure  itself  may  range  from  cracks,  detectable  only  by  microscopic 
examination,  to  sublimation  or  disintegration  of  part  of  the  material.  The  latter 
type  of  failure  is  referred  to  as  ultimate  or  cohesive  spall;  the  former  type 
represents  a  lower  limit  to  a  different  and  less  drastic  failure  mode  called 
ductile  spall.  In  addition,  materials  which  exhibit  a  pressure-induced  phase 
transformation  may  fail  in  an  intermediate  mode  known  as  phase  transformation 
spall.  Since  both  ultimate  and  phase  transformation  spall  occur  at  relatively 
high  stress  levels  compared  to  ductile  spall,  and  since  they  have  received  only 
a  small  amount  of  experimental  or  theoretical  attention,  they  will  be  mentioned 
only  briefly  in  this  study.  The  major  emphasis  will  be  placed  on  the  more 
practically  interesting  failure  mechanism  known  as  ductile  spall. 
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Ultimate  or  Cohesive  Spall.  This  mode  of  spall  has  been  obierved  by  a 
few  investigators,  most  notably  McQueen  and  Marsh. ^  By  generating  ihort  dura¬ 
tion  pulses  with  an  explosively-driven  flyer  plate,  and  assuming  hydridynamic 
material  response,  the  ultimate  strength  of  single-crystal  and  polycr>atalline 
copper  was  found  to  be  in  excess  of  150  kbar.  The  distinguishing  feature  of  the 
observed  failure  mode  was  that  the  material  literally  disintegrated  or  sublimed, 
in  contrast  to  the  behavior  at  lower  spall  thresholds  where  flaw-related  mechanisms 
such  as  crack  nucleation  and  propagation  are  prominent.  It  is  interesting  to  note 
the  rather  favorable  comparison  between  the  experimentally  obtained  ultimate  spall 
threshold  and  the  static  theoretical  strength  of  the  material  derived  from  con¬ 
sideration  of  the  interatomic  forces  in  a  crystal.  Using  the  equation  for  the 
theoretical,  stati:  cohesive  of  .05E  -  .10E,  and  assuming  a  modulus  of  18  x  10**  psl 
for  copper,  the  strength  is  found  to  lie  between  62  and  125  kbar. 

Phase -Transformation  Spall.  This  mode  of  failure  occurs  only  in  materials 
which  exhibit  a  pressure-induced  phase  transformation.  It  has  been  observed  in 
Armco  iron  by  Erkman^  and  Moss  and  Glass.''  '  Both  of  these  investigations  used 
plate-slap  impact  techniques  to  induce  stresses  in  excess  of  the  phase “transforma¬ 
tion  pressure  (130  kbar),  thereby  creating  a  characteristic  double  plastic  shock 
wave  structure.  The  phase- transformation  spall  threshold,  which  is  the  negative 
of  the  phase- transformation  pressure,  does  not  represent  a  point  of  change  in 
failure  mechanism,  but  rather  marks  a  point  of  change  in  the  appearance  of  the 
spall  surface.  Thi3  change  is  attributed  to  the  phase  transformation.  As  induced 
pressures  are  increased  through  the  phase- transformation  pressure,  the  spall 
plane  takes  on  a  characteristically  smooth  appearance,  although  the  failure 
mechanism  remains  one  of  crack  nucleation  and  propagation.  Phase-transformation 
spall  has  also  been  observed  in  steel  by  Novikov  and  others, and  in  iron  and 
steel  by  Ivanov  and  others both  studies  used  contact  high  explosives  to 
induce  the  required  pressures. 

Ductile  Spall.  Low-pressure  or  ductile  spall  involves  the  same  basic 
mechanisms,  i.e.,  crack  nucleation  and  propagation,  that  are  associated  with  lower 
strain-rate  failure.  Three  separate  levels  of  ductile  spall  have  been  recognized 
and  classified  according  to  the  severity  of  damage.  In  order  of  increasing  damage 
they  are;  incipient,  intermediate  >.r.d  complete. 
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The  incipient  spall  threshold  is  defined  as  that  combination  of  stress 
and  time  (pulse  duration)  below  which  no  damage  to  the  specimen  would  be  visible. 

In  order  to  determine  the  visibility  of  failure,  it  is  necessary  to  section, 
polish,  and  etch  the  specimen  and  perform  a  metallographic  examination  at  about 
100  x  magnification.  The  large  number  of  parameters  influencing  the  spall  thresh¬ 
old  precludes  a  general  description  of  the  appearance  of  the  fractured  region, 
although  it  can  be  stated  that  several  usually  transgranular  microcracks  or  voids 
with  dimensions  on  the  order  of  a  grain  diameter  would  be  evident.  Also,  the 
cracks  would  be  generally  parallel  and  would  lie  in  a  fairly  narrow  hand  parallel 
to  the  free  surface  of  the  specimen.  Photomicrographs  illustrating  the  appearance 
of  incipient  spall  in  6061-T6  aluminum  and  copper  are  given  in  References  12  and 
13,  respectively. 

If  the  tensile  stress  and/or  puls«  duration  is  greater  than  that  required 

for  incipient  spall,  the  degree  of  damage  is  correspondingly  greater;  i.e.,  the 

cracks  are  more  numerous  and  larger,  and  some  are  joined  together.  A  criterion 

for  quantitative  determination  of  the  degree  of  damage  has  been  suggested  by 

Herrmann  of  Sandia,  who  proposed  that  the  damaged  specimens,  including  the  spall 

plane,  be  cut  into  tensile  members  and  pulled  apart  in  a  tension  testing 
(13) 

machine.  In  this  way,  a  quantitative  measure  of  the  degree  of  damage  is  pro¬ 

vided  by  the  residual  strength  of  the  tensile  member.  As  the  tensile  stress 
and/or  pulse  duration  is  increased  beyond  the  incipient  spall  threshold  values, 
the  degree  of  damage  to  the  specimen  increases  and,  hence,  its  residual  strength 
decreases.  At  some  point,  however,  the  residual  strength  of  the  specimens  ceases 
to  decrease  with  increase  in  the  severity  of  loading.  The  tensile  stress  and 
pulse  duration  corresponding  to  this  point  are  said  to  define  the  intermediate 
spall  threshold. 

It  a  specimen  is  subjected  to  a  tensile  pulse  of  sufficient  magnitude 
and/or  duration  in  excess  of  the  intermediate  spall  threshold  values,  the  damage 
will  be  so  great  that  a  fairly  intact  piece  of  material  will  separate  or  be 
detached  from  the  specimen.  The  combination  of  stress  and  pulse  duratic  .  at  which 
this  occurs  is  known  as  the  complete  spall  threshold.  This  threshold  is  of  greater 
practical  importance  thar  the  others  and  has,  historically,  been  the  object  of 
che  greatest  study,  e.g.,  Hopkinson,  Kolsky,  Rinehart.  Obviously,  the  residual 
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strength  of  the  material  above  this  threshold  is  zero,  and  a  complete  charac¬ 
terization  of  the  three  ductile  spall  thresholds  in  terms  of  residual  strength 
would  appear  as  shown  in  Figure  7 . 


FIGURE  7.  RESIDUAL  STRENGTH  OF  TARGETS  USING  1.5  Ps  SQUARE  PULSES 

*  The  24-kbar  data  were  obtained  with  a  grooved  target. 

Edge  effects  prevented  complete  separation  at  22  kbar. 

Photomicrographs  illustrating  the  appearance  of  the  fracture  zone  in  copper  at 
the  three  ductile  spall  thresholds  are  given  ia  Reference  13. 

Some  of  the  more  recent  investigators  of  the  complete  spall  threshold 
have  used  specimen  configurations  specifically  designed  to  eliminate  the  influence 
of  edge  effects  and  shear  failure,  which  necessarily  enter  into  the  complete  spall 
of  a  standard  plate  specimen-  By  using  tapered  plugs  fitted  into  plates  of 
similar  material,  Smith obtained  complete  spall  thresholds  for  copper  which 
were  lower  than  those  obtained  with  standard  specimens,  tfhen  a  scheme  such  as 
this  is  employed,  however,  care  must  he  taken  to  ensure  that  uniaxial  strain 
conditions  are  maintained  while  the  material  i3  in  tension.  If  this  is  not  done, 
the  different  stress  states  occurring  in  compression  and  tension  must  be  accounted 
for  in  the  solution  of  the  wave  propagation  equations.  Glass  and  others, 
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experimenting  with  single  crystals  of  copper,  used  cylindrical  specimens  of  L/D  >  1 
embedded  in  copper  plates.  They  observed  rather  severe  necking  in  the  region  of 
the  specimen  adjacent  to  the  spall  zone,  indicating  a  stress  state  similar  to  that 
prevailing  in  a  standard  tensile  test,  i.e.,  uniaxial  stress. 

Parameters  Influencing  Ductile  Spall  Thresholds 

A  sufficient  amount  of  experimental  and  theoretical  evidence  has  now  been 
gathered  to  support  the  assertion  that  the  "critical  normal  fracture  stress" 
criterion  for  spall  must  be  at  lease  extended,  if  not  totally  revised.  This 
criterion,  which  states  that  a  material  will  spall  instantaneously  and  completely 
when  a  unique  critical  value  of  normal  tensile  stress  is  attained,  has  yielded 
erroneous  and  even  contradictory  predictions  of  spall  in  several  instances.  Indica¬ 
tions  are  that  a  realistic  spall  criterion  must  include  at  least  one,  if  not  several 
more  parameters  to  account  for  the  anomalies  in  experimentally  observed  spall 
behavior.  In  this  section,  those  parameters  which  have  been  found  to  affect  spall 
behavior  will  be  enumerated  and  discussed,  with  reference  to  specific  experiments. 

Time-Related  Parameters.  One  of  the  first  indications  that  a  critical 
normal  fracture  stress  did  not  constitute  an  adequate  spall  criterion  was  the 
disparity,  both  in  numbers  and  sizes,  between  predicted  and  measured  scabs  in 
multiple-spall  experiments.  In  his  experiments  with  several  metals  and  glass 
subjected  to  explosively  induced  pressures,  Broberg^"^  attributed  these  discrepancies 
to  the  nonins tantaneous  fracture  behavior  of  ductile  materials.  If  a  material  which 
is  subjected  to  a  triangular  compressive  pulse  does  not  spall  instantaneously  when 
the  net  reflected  tension  reaches  a  critical  value,  then  the  magnitude  of  the  pulse 
available  to  cause  additional  spalling  upon  reflection  from  the  first  spalled 
surface  is  reduced. 

In  a  later  article, reporting  on  experiments  done  with  granite  and 
Incite,  Broberg  explained  the  spall  delay  time  in  terms  of  the  mechanisms  associated 
with  fracture  at  lower  strain  rates,  i.e.,  crack  nucleation  and  propagation.  Thus, 
the  stress-dependent  fracture  delay  time  decreases  with  increasing  stress,  not 
only  because  more  cracks  are  nucleated,  but  also  because  the  acceleration  of  the 
crack  tips  increases.  Similarly,  ductile  materials  require  more  time  to  fracture 
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than  brittle  materials  because  stresses  are  relieved  not  only  by  fracture,  but 
also  by  flow.  Ductility  limits  the  magnitude  of  stress  concentrations  and  causes 
the  tips  of  cracks  to  become  blunt,  thereby  retarding  both  the  velocity  and  the 
acceleration  of  crack  propagation. 

Results  of  plate-slap  tests  provide  additional  evidence  for  the  necessity 
of  including  a  time-related  parameter  in  the  spall  criterion.  Butcher  and  others 
point  out  that,  for  the  same  material,  the  difference  between  the  spall  threshold 
stress  generated  by  contact  explosives  (triangular  pulses)  and  the  spall  threshold 
stress  obtained  by  plate-slap  impact  (rectangular  pulses)  indicates  a  time  dependence 
for  spall  failure.  In  the  cases  where  a  strict  comparison  was  possible,  the 
explosively  induced  spall  threshold  stress  was  always  greater  than  the  plate-slap 
spall  threshold. 

(181 

In  plate-slap  experiments  on  6061-T6  aluminum,  Blincow  and  Keller'  ' 
observed  that  the  flyer  velocity  required  to  cause  incipient  spall  decreased  with 
increasing  flyer  thickness  in  an  approximately  linear  manner.  Since  impact 
velocity  can  be  directly  related  to  stress,  and  flyer  thickness  can  be  similarly 
related  to  pulse  duration,  a  definite  time  dependence  of  spall  is  implied.  Another 
interesting  observation  made  by  these  investigators  was  that  the  thickness  or  extent 
of  the  spalled  region  increased  with  decreasing  impact  velocity  (above  the  incipient 
spall  threshold).  This  phenomenon  can  also  be  attributed  to  a  time-dependent  spall 
mechanism  by  recognizing  that  for  a  given  pulse  width,  if  the  stress  is  greater 
than  that  required  for  spall,  the  material  will  fail  more  quickly  than  it  would  at 
the  threshold  value  of  stress.  Due  to  this  more  rapid  spall  at  higher  stress  levels, 
the  remainder  of  the  tension  wave  reflects  sooner  from  the  spalled  surface  as  a 
compressive  wave  to  prevent  additional  fracturing  near  the  original  spall.  Also, 
since  the  width  of  the  pulse  that  is  transmitted  across  the  spall  plane  prior  to 
spall  is  less  than  that  at  lower  stress  levels,  the  material  between  the  impact 
surface  and  the  spall  plane  is  less  likely  to  experience  subsequent  fracture 
(multiple  spall). 

Other  investigators  have  found  evidence  for  time- dependent  spall  in 

(13) 

experiments  with  copper.  In  his  plate-slap  experiments,  Smith'  found  that  the 
required  tensile  stress  decreased  with  increasing  pulse  duration  for  all  ductile 
spall  thresholds.  Of  additional  interest  in  this  reference  are  several  photo¬ 
micrographs  illustrating  spall  damage  for  different  combinations  of  stress  and  pulse 
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duration.  These  pictures  show  increasing  damage,  i.e.,  larger  cracks  and  voids, 
for  the  case  when  stress  is  held  constant  and  pulse  duration  increased  and  the 
case  when  pulse  duration  is  kept  constant  and  stress  is  increased. 

(14) 

Glass  and  others  used  contact  explosives  to  study  the  fracture  of 
copper  single-crystals.  These  authors  point  to  the  presence  of  a  necked-down 
region  adjacent  to  the  spall  zone  as  definite  proof  that  the  material  does  not 
fracture  instantaneously. 

Not  all  investigators  explicitly  introduce  time  as  a  parameter  with 

(19) 

which  to  augment  a  spall  criterion.  Al'Tshuler  and  othersv  '  attribute  the 
difference  between  explosive  and  plate-slap  spall  threshold  stresses  to  the 
different  tensile  strain  rates  corresponding  to  these  loadings.  It  is  interesting 
to  note  that,  contrary  to  the  results  of  the  comparison  made  by  Butcher, (17)  they 
found  the  spall  threshold  stress  for  copper  obtained  by  plate-slap  impact  to  be 
more  than  twice  the  explosively  induced  threshold  stress.  They  estimate  the  strain 
rates  as  105  sec  1  for  explosive  loading  and  107  sec  1  for  plate  slap. 

Whiteman has  introduced  tensile-stress  unloading  rate  to  account  for 
the  nonunique  spall  threshold  stresses  which  he  measured  in  aluminum,  mild  steel, 
brass,  and  copper  by  plate-slap  impact  techniques.  Plausibility  of  the  results  is 
indicated  to  the  extent  that  curves  of  spall  stress  versus  this  stress  rate 
generally  show  spall  stress  approaching  the  theoretical  strength  at  the  highest 
stress  rates.  In  addition,  the  trends  of  the  curves  appear  to  be  consistent  with 
theoretical  considerations.  If  it  is  assumed  that  cracks  are  nucleated  by  dis¬ 
location  pile-ups  at  grain  boundaries,  then  the  stress  necessary  to  initiate 
fracture  must  increase  with  increasing  stress  rate  (and,  hence,  shorter  load- 
application  times)  to  give  the  dislocations  sufficient  velocity  to  reach  the 
boundary.  More  discussion  of  the  stress-rate  dependent  spall  criterion  will  be 
given  in  the  section  devoted  exclusively  to  spall  models  (pages  23-32). 

Other  Parameters.  In  contrast  to  time-related  effects,  evidence  per¬ 
taining  to  the  influence  of  other  parameters  on  ductile  spall  thresholds  is 
comparatively  meager.  Furthermore,  some  of  this  evidence  is  not  only  rather 
tenuous,  but  contradictory  as  well.  For  the  sake  of  completeness,  however,  results 
pertaining  to  the  influence  of  all  parameters  (other  than  time,  of  course)  will  be 
presented  here. 
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The  previously  discussed  discrepancies  observed  in  multiple-spall 

experiments  evoked  an  explanation  other  than  that  of  a  time- dependent  spall 

(21) 

mechanism.  Rinehart  felt  that  the  discrepancies  could  be  explained  by  the 

fact  that  not  all  spall  planes  experienced  the  same  peak  compressive  stress  before 

going  into  tension.  In  other  words,  the  critical  normal  fracture  stress  of  a 

material  depends  on  the  transient  compressive  stress  immediately  preceding 

tension.  To  the  extent  that  explosively  generated  pulses  can  be  controlled,  and 

their  shapes  accurately  described,  Rinehart's  suspicions  are  borne  out  by  the 

(22) 

experiments  of  Buchanan  and  James  with  mild  steel.  Varying  the  magnitude  of 
the  precompression  pulse  by  using  specimens  of  different  thicknesses,  they  found 
that  the  critical  normal  fracture  stress  increased  in  an  approximately  linear 
manner  with  increasing  precompression.  In  particular,  the  tensile  stress  required 
to  cause  rpall  increased  by  about  SO  percent  when  the  magnitude  (peak  stress)  of 

the  preceding  compressive  pulse  was  tripled.  However,  n  these  experiments,  the 
rate  of  stress  application  also  increased  with  increasing  peak  stress  so  that  the 
precompression  explanation  of  the  observed  effect  is  not  unique. 

More  work  should  be  done  to  determine  what  effects  the  compression  phase 
of  the  stress-time  history  has  on  the  spall  behavior  of  materials,  if  for  no  other 
reason  than  the  fact  that  crack  nucleation,  which  is  dependent  on  shear  stresses 
rather  than  normal  stresses,  can  occur  in  the  compressive  as  well  as  the  tensile 
phase.  In  both  phases,  the  stress  states  are  nonhydrostatic.  Note  that  the  trend 
of  the  results  described  above,  i.e.,  increasing  strength  with  increasing  pre¬ 
compression,  runs  contrary  to  what  would  be  expected  if  Increased  precompression 
implied  a  greater  amount  of  crack  nucleation. 

The  influence  of  ambient  temperature  on  spall  thresholds  is  another  area 

warranting  further  investigation.  For  the  rather  extensive  range  of  temperatures 

(23^ 

they  considered,  O'Brien  and  Davis v  '  found  no  change  in  the  fracture  stress  of 

high-purity  polycr^  italline  aluminum  subjected  to  explosive  loads.  In  contrast  to 

(12) 

this  are  the  results  obtained  by  Penning  and  others  from  plate-slap  tests  on 
lucite.  With  a  very  limited  number  of  tests,  covering  only  a  small  temperature 
range,  they  found  that  the  fracture  stress  increased  with  increasing  temperature. 

It  appears  that  neither  of  the  above-mentioned  investigators  accounted 
for  the  influence  of  temperature  on  the  equation-of-state,  and  consequently,  on 


1 
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the  wave  propagation.  Obviously,  this  type  of  effort  must  accompany,  if  not  pre¬ 
cede,  any  investigation  into  the  effect  of  ambient  temperature  on  spall  thresholds. 

(241 

The  results  of  a  study  done  by  Places!  and  Watt  show  that  yield 
strength  can  play  a  double  role  in  the  spallation  phenomenon.  In  addition  to 
affecting  the  shape  of  the  stress  pulse  and,  therefore,  the  stress-time  history 
at  the  spall  plane,  the  yield  strength  of  a  material  influences  the  amount  of 
bulge  at  the  rear  surface  of  a  specimen  that  has  not  completely  spalled.  Thus, 
with  specimens  not  designed  to  be  free  of  edge  effects  and  shear  failure,  such  as 
described  previously,  yield  strength  would  influence  both  the  intermediate  and 
complete  spall  thresholds. 

In  this  study,  the  effect  of  yield  strength  on  rear  surface  bulge  was 

demonstrated  experimentally  by  impacting  7075-T6  aluminum  targets  at  various 

temperatures  (and,  consequently,  different  yield  strengths)  with  1/4-inch-diameter 

aluminum  spheres.  The  bulge  distance  versus  impact  velocity  curves  for  7075-T6 

aluminum  at  500°  F  and  7075-0  aluminum  at  72°  F,  which  have  equal  yield  strengths, 

were  approximately  identical.  This  result  indicated  that  temperature  did  not 

(23) 

affect  spall,  in  agreement  with  the  conclusions  of  O'Brien'  7  discussed  previously. 
Thus,  the  results  obtained  in  tests  at  lower  temperatures,  i.e.,  smaller  bulge 
distances  for  all  impact  velocities,  could  justifiably  be  attributed  to  the 
influence  of  yield  strength  alone. 

The  necessity  for  investigating  the  influence  on  spallation  of  micro¬ 
scopic  material  properties,  such  as  structure  and  composition,  has  been  recognized, 
and  some  work  has  been  done.  In  the  study  previously  mentioned, (23)  O'Brien  and 
Davis  found  that  the  fracture  stress  remained  constant  for  aluminum  single  crystals 
at  various  orientations,  zone-refined  polycrystalline  aluminum,  polycrystalline 
aluminum  rods  that  had  been  cold  worked  to  various  textures,  and  2024  aluminum- 
copper  alloy  in  the  hardened  and  solution-treated  condition.  Although  the  ap¬ 
pearance  of  the  spall  was  influenced  by  the  ductility  of  the  specimens,  the 
conclusion  reached  was  that  the  structural  defect  causing  failure  was  present  in 
an  as-grown  crystal  or  was  caused  by  the  precompression  shock. 

In  their  experiments  with  lead,  aluminum,  and  copper,  subjected  to 

(25) 

explosive  loads,  Vitman  and  others  found  spall  threshold  stresses  to  be 
insensitive  to  alloying. 
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The  structural  parameter  that  nas  been  found  to  have  the  greatest 
Influence  on  spallation  is  the  material  form  (I.e.,  bar  or  plate).  Using  plate- 
slap  techniques  to  study  spallation  in  copper,  Plauson  and  others '  '  found  that 

the  spall  stresses  for  annealed  bar  s;  ’imens  (grains  aligned  normal  to  the  speci¬ 
men  faces)  were  more  than  twice  the  stresses  for  plate  specimens  (grains  aligned 
parallel  to  the  specimen  faces).  In  addition,  the  spall  stresses  of  the  plate 
specimens  were  j.ound  to  be  relatively  insensitive  to  material  condition  (i.e., 
annealed,  1/4  hard,  or  hard). 

The  individual  cracks  constituting  the  entire  spall  zone  were  observed 
to  be  of  a  t r an s granular  nature  in  both  the  bar  and  plate  specimens,  but  the  gross 
appearance  of  the  spall  differed  greatly.  In  the  plate  material,  the  individual 
fractures  were  cleavage  fractures  which  showed  little  attendant  grain  deformation, 
while  in  the  bar  material,  shear  fractures  and  large  amounts  of  grain  deformation 
were  observed.  In  both  bars  and  plates,  the  gross  spall  formed  on  a  plane  per- 
pendicular  to  the  direction  of  impact.  However,  the  individual  cracks  were  oriented 
in  the  direction  of  grain  alignment. 

Similar,  but  less  quantitative  results  on  the  effect  of  grain  alignment 

(13) 

in  copper  have  been  obtained  by  Smith.  This  reference  contains  excellent 

photomicrographs  which  illustrate  the  marked  differences  in  the  appearance  of  the 
spall  zones  in  bar  and  plate  specimens. 

In  a  series  of  plate-slap  experiments  on  6061-T6  aluminum  plate  stock, 
Butcher^)  found  that  the  critical  values  of  stress  and  pulse  duration  required 
for  crack  nucleation  were  insensitive  to  the  angle  between  the  impact  velocity 
vector  and  the  plane  of  grain  elongation  and  impurity  stratification.  Crack 
propagation,  on  the  other  hand,  was  found  to  be  extremely  sensitive  to  this 
orientation.  Cracks  usually  propagated  in  the  plane  of  weakness,  i.e.,  grain 
elongation  and  impurity  stratification,  regardless  of  the  plane's  orientation. 

Thus,  when  this  plane  was  parallel  to  the  direction  of  impact,  little  propagation 
was  noticed  because  the  tension  normal  to  this  plane  is  neither  uniform  nor  of 
maximum  magnitude.  In  this  crientation,  evidence  of  massive  plastic  flow  was 
found  in  the  area  of  the  spall  zone. 
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It  appears  that  the  increased  spall  strength  of  specimens  in  which 
grain  elongation  is  parallel  to  the  direction  of  impact  is  due  to  the  fact  that 
cracks  nucleated  in  these  specimens  are  not  aligned  perpendicular  to  the  direction 
of  maximum  tension  and,  thus,  do  not  propagate  as  readily  nor  coalesce  as  easily 
as  cracks  aligned  parallel  to  the  spall  plane.  Since  incipient  spall  involves 
verv  '  '  —ack  propagation,  the  insensitivity  of  this  threshold  to  grain  align- 

/  l  o\ 

ment.  ”  /ed  by  Butcher,  is  not  surprising.  Blincow  and  Keller  '  1  also 

repoi  e  incipient  spall  threshold  of  6061-T6  aluminum  is  the  same  for 

both  plu  -nd  bar  specimens,  but  certainly  more  data  is  needed. 

These  observations,  and  the  remarks  made  previously  in  connection  with 
the  effect  of  precompression,  indicate  that  significant  progress  toward  under¬ 
standing  the  spallation  phenomenon  must  necessarily  involve  a  deeper  under¬ 
standing  of  the  roles  played  by  the  basic  fracture  mechanisms. 

Spall  Models 

Except  for  very  simple  pulse  shapes,  analytical  predictions  of  spall 
(i.e.,  thickness  and  number)  are  handled  mo3t  efficiently  by  a  computer  in 
conjunction  with  wave  propagation  calculations.  The  critical  normal  fracture 
stress  spall  criterion,  being  the  simplest  and  requiring  the  least  data,  is 
presently  being  used  in  several  codes  and  requires  no  further  discussion.  The 
more  sophisticated  models,  incorporating  such  parameters  as  time,  rate,  tempera¬ 
ture,  and  thus  requiring  more  experimental  data,  are  at  various  stages  of  develop¬ 
ment.  Some  have  had  successful,  if  limited,  application  and  verification,  while 
others  remain  at  the  purely  conceptual  stage.  With  one  exception,  which  will  be 
discussed  separately,  all  of  these  spall  models  are  phenomenological  in  nature; 
e.g.,  spall  stress  is  related  to  time,  rate,  etc.,  by  experimentally  determined 
constants  which  characterize  the  spall  behavior  of  the  material.  The  explanatory 
model,  based  on  the  principle  of  bond  rupture,  represents  an  attempt  to  predict 
spall  behavior  from  consideration  of  the  physical  mechanisms  actually  involved. 

Phenomenological  Spall  Models.  Aside  from  being  a  numerical  application 
of  a  particular  spall  criterion,  a  spall  model  should  also  have  the  capability  to 
predict  the  spall  behavior  of  a  material  subjected  to  c  pulse  of  arbitrary  shape. 
Thus,  if  an  experimentally  determined  relationship  between  fracture  delay  time 
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and  consfant  applied  tensile  stress  is  known,  the  spall  model  must  contain  some 

statement  as  to  how  this  information  is  to  be  used  when  the  material  is  subjected 

(121 

to  a  variable  stress.  Penning  and  others  suggest  that  the  spall  delay  time 
of  a  material  under  variable  tensile  stress  is  cumulative,  i.e.,  the  material  will 
spall  at  that  time,  T,  given  by 


rT  dt 
j  r 


i  ; 


o  s 

where  t  ■  t  (a)  is  the  delay  time  versus  constant 
8  8 

tensile  stress  relationship. 

To  see  this  more  clearly,  consider  the  sketch  of  the  stress-time  history,  a 
shown  in  Figure  8. 


(3) 


-  o  (t). 


FIGURE  8.  ASSUMED  TENSILE  STRESS  TIME-HISTORY 
AT  THE  SPALL  PLANE 


At  time  t  ■  t,  spall  occurs,  i.e., 
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It  can  be  seen  thAt  spall  does  not  jccur  until  the  percent  of  damage  associated 
with  the  individual  approximating  pulses  (aB,  At„)  accumulates  to  100  percent. 

The  obvious  shortcoming  of  this  scheme  is  that  the  effect  of  previous 

stressing  is  r.ot  included  in  the  t  versus  a  relationship.  For  each  time  incre- 

8 

ment,  Atn,  the  corresponding  delay  time,  t  (cB),  is  based  on  an  being  applied 

s 

instantaneously,  with  no  previous  stress  present.  From  these  considerations,  it 
would  appear  that  the  t  obtained  by  this  method  wculd  be  too  large. 

Butcher  and  others  have  utilized  the  cumulative  method  to  obtain 
closed  fora  expressions  y.r  T.  By  considering  only  triangular  pulses,  and  as¬ 
suming  a  definite  functional  relationship  between  t  and  a,  they  were  able  to 

s 

integrate  Equation  3  and  solve  for  t  in  terms  o*  a  material  spatial  coordinate, 
the  shock  velocity,  and  parameters  characterizing  the  pulse  shape  and  the  t 

8 

versus  o  relationship. 

/OO  \ 

Tuler  and  Butcher'  '  incorporated  the  cumulative  spall  model  into  an 

elastic-plastic-vavc  propagation  code.  They  used  a  t  versus  a  relationship  of 

8 

the  fora 

t  -  k(a-a)"X  ;  (5) 

s  o 

where  k  and  X  are  positive  constants  and 
ao  represents  the  minimum  tensile  stress 
required  for  failure,  regardless  of  time 
duration. 

Predictions  made  with  this  model  were  compared  to  spall  data  for  1100  aluminum 
which  had  been  subjected  to  explosively  induced  pulses. 

More  fundamental  difficulties  are  encountered  when  the  spall  delay  time 
versus  applied  tensile  stress  relationship  is  examined.  Two  of  the  more  common 
forms  of  this  relationship  are  given  by 


t  -  Xe"^ 
s 

(References  29  and  30), 

(6) 

t  *  acP 

s 

(Reference  31); 

(7) 

where  X,  p, 

a,  and  @  are  constants. 

Equation  6  represents  not  only  the  experimentally  observed  behavior  of 
a  wide  variety  of  materials  when  tg  is  relatively  large  (milliseconds  and  above), 
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but  has  been  derived  by  Waldorf from  theoretical  considerations  and  shown 
to  hold  for  al’  values  of  tg.  In  this  theory,  which  will  be  discussed  sub¬ 
sequently  in  more  detail  (see  page  28),  the  constants  X  and  p>  are  identified 
with  such  physical  parameters  as  sublimation  energy,  Boltzmann's  constant, 
absolute  temperature  Inside  the  shock,  and  the  lattice  vibrational  frequency. 

Equation  7  merely  represents  a  fit  to  data  gathered  for  several 
materials,  most  of  which  were  ablators  and  plastics. 

Regardless  of  the  fora  of  the  delay  time  versus  constant  applied-stress 
relationship,  the  difficulty  of  applying  the  relationship  to  the  prediction  of 
spall  when  a  material  is  subjected  to  unsteady  tensile  stress  is  one  of  logical  con¬ 
tradiction.  When  the  stress  causing  spall  is  not  large  compared  to  the  dynamic 
yield  stress,  the  hydrodynamic  interpretation  of  material  response  is  not 
justified.  Elastic-plastic  wave-propagation  calculations  show  that  the  stress¬ 
time  history  at  the  spall  plane  is  very  unsteady  due  to  the  complicated 
interaction  of  the  unloading  compressive  wave  and  the  elastic-plastic  rare¬ 
faction  from  the  free  surface.  Thus,  data  describing  a  spall  threshold  in  term- 
of  stress  and  pulse  duration  necessarily  imply  that  the  fine  structure  of  the 
actual  stress-time  history  at  the  spall  plane  is  unimportant,  and  that  spall 
be  characterized  in  terms  of  an  equivalent  rectangular  pulse.  This  equivalen: 

rectangular  pulse  can  be  obtained  in  various  ways.  In  his  experiments  with 
(13) 

copper,  Smithv  completely  neglected  the  effect  of  the  elastic  precursor. 

(12) 

Others  have  performed  sophisticated  calculations  of  the  stress-time  history 
at  the  spall  plane,  and  then  they  appear  to  have  "eyeballed"  an  equivalent 
rectangular  pulse  with  which  to  describe  the  spall  threshold. 

From  these  constat- ations ,  it  must  be  concluded  that  the  only  con¬ 
sistent  way  to  predict  the  spall  behavior  of  a  material  subjected  to  an  arbitrary 
finely-structured  pulse,  from  this  type  of  data,  is  to  idealize  the  actual  pulse 
in  a  manner  identical  to  that  which  was  used  to  obtain  the  data.  From  a  physical 
standpoint,  a  more  desirable  approach  would  be  to  introduce  additional  parameters 
to  characterize  the  actual  shape  of  the  spall-threshold  pulse.  Again,  since 
fracture  may  be  considered  basically  a  mechanism  of  energy  transfer,  it  may  be 
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more  realistic  to  describe  spall  in  terms  of  total  mechanical  energy  introduced 
to  the  spall  plane.  A  spall  criterion  based  on  this  concept  would  have  the  form 

k  «  Jovdt  ;  (8) 

o 

where 

k  =  an  empirically  determined  constant  in  units  of  energy  per 
unit  area, 

t  «  time  at  which  spall  occurs , 

a  m  instantaneous  net  tensile  stress  at  the  spall  plane, 
v  *  instantaneous  net  particle  velocity  at  the  spall  plane, 
t  ■  the  length  of  time  for  which  the  spall  plane  has  been 
in  tension. 

If  only  rectangular  pulses  are  considered,  and  if  linear  elastic  material  approxi¬ 
mations  are  made,  i.e.,  a  *  pcv,  then  this  equation  can  be  integrated  to  yield 

t  -  t  -  kpco‘2  .  (9) 

s 

Note  the  formal  similarity  to  Equation  7  which  represents  a  fit  to  experimental 
data.  Equation  7  also  involves  a  negative  exponent,  although  it  is  of  different 
magnitude. 

A  spall  model  based  on  the  criterion  of  stress-rate -dependent  spall 

(22) 

stress  has  been  used  successfully  by  Breed  and  others  at  Loo  Alamos  Scientific 

f2Gv 

Laboratories  (LASL) .  These  investigators  found  that  Whiteman's v  7  data,  relating 
apall  stress  to  the  tensile  stress  unloading  rate  for  aluminum  and  copper,  could 
be  put  in  the  computationally  convenient  form  of 


where 


(10) 


ag  *  spall  stress, 

Arr 

.gX  *■  instantaneous  net  tensile-stress  unloading  gradient  at  spall  plane, 

A,B  =  experimentally  determined  constants. 

Using  this  model  in  the  LASL  hydrodynamic  SIN  code,  the  damage  done  to 
a  2.5  cm  thick  aluminum  block,  subjected  to  the  contact  explosion  of  a  10  cm 
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block  of  composition  B-3,  was  predicted  to  be  10  spalls  in  23  p.8.  A  radiograph 
of  an  actual  test  specimen  taken  at  25  Us  showed  II  ^  1  spalled  layers.  Close 
agr“°ment  between  predictions  of  the  code  and  experimental  results  was  also 
found  for  the  case  of  a  copper-aluminum-copper  sandwich  subjected  to  plate-slap 
impact. 

These  results  illustrate  how  successful  a  phenomenological  spall  model 
can  be,  provided  that  consistency  is  maintained  between  the  analysis  of  test  data 
and  analytical  predictions  of  material  behavior  based  on  this  data.  Since  it  is 
now  generally  accepted  that  the  hydrodynamic  material  approximation  yields 
erroneous  results,  except  at  extremely  high  pressures,  a  more  physically  satis¬ 
fying  approach  would  be  to  incorporate  this  spall  model  into  a  code  for  elastic- 
plastic  stress  wave  calculations. 

Explanatory  Spall  Model.  The  spall  model  developed  by  Waldorf^"^  is 

based  on  the  fundamental  physical  criterion  that  atomic  bonds  must  be  broken  in 

order  for  fracture  to  occur.  By  using  this  criterion  in  conjunction  with  several 

material  idealizations,  Waldorf  was  able  to  derive  a  relationship  between  spall 

delay  time,  t  ,  and  applied  tensile  stress,  a  ;  for  homogeneous  materials,  i.a., 
s  o 

single  crystals.  This  relationship  is  then  extended  and  modified  to  describe  the 
fracture  of  materials  in  which  crack  propagation  plays  an  important  role,  i.e., 
polycrystalline  materials.  The  essential  features  involved  in  the  development 
of  these  relationships  are  outlined  in  this  section. 

For  a  microscopically  homogeneous,  isotropic,  and  flawless  material 
loaded  in  uniaxial  strain  with  the  bonds  parallel  to  the  applied  stress  carrying 
the  entire  load,  the  rate  at  which  bonds  are  broken  can  be  written  as 


dn 

dt 


SSL  *  -n?v 


(11) 


where 


n  *  the  number  of  unbroken  bonds  per  unit  area  perpendicular  to 
the  stress, 

P  *  the  probability  of  an  atomic  bond  rupturing  per  unit  time, 
b 

(This  is  a  function  of  the  applied  stress  and  temperature.) 


£  ' 


i  | 

.. 


-»  i 
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Now  assuming  that  the  atoms  move  about  their  equilibrium  positions  in  the  lattice 
with  a  Maxwellian  velocity  distribution  in  the  direction  of  the  applied  stress, 
and  defining  bend  rupture  as  occurring  when  an  atom  has  traveled  a  distance  equal 
to  the  equilibrium  spacing,  6q,  the  rupture  probability,  P^,  is  derived  as 

Pb  -  iue"E/kt  ;  (12) 

where 

a)  «  the  vibrational  frequency  of  the  lattice, 
k  «*  Boltzmann's  constant, 

T  *  absolute  temperature, 

E  *  the  energy  per  bond  with  which  the  atoms  are  bound  into 
the  lattice. 

In  the  unstressed  state,  the  binding  energy,  E,  is  equal  to  Eq,  which  is  approxi¬ 
mately  the  sublimation  energy. 


The  dependence  of  E  on  applied  stress,  is  found  by  assuming  that  E 
varies  linearly  with  lattice  spacing,  6,  and  by  using  an  average  value  of  bulk 
modulus  to  relate  a0  to  6.  By  further  assuming  that  fracture  occurs  at  a  time 
equal  to  the  mean  lifetime  of  the  bonds  and  that,  during  this  time,  the  number  of 
unbroken  bonds  may  be  averaged  to  a  constant  value,  the  E  versus  Oo  relationship 
becomes 


where 


i 
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When  oQ  is  constant,  this  reduces  to 

*,  ■  « 


6  a 


KT 


(16) 


A  less  approximate  approach  would  be  to  relate  E  and  6  by  the  Morse 
potential  function,  i.e., 


where  a  is  a  constant  such  that  E  ■  o  when  6  »  26Q,  i.e.,  the  atoms  become  unbound 
when  forced  apart  by  an  amount  6Q.  Then,  using  an  average  value  of  effective 
modulus  in  uniaxial  strain,  E,  it  can  be  written  that 


and 


(18) 


E  “  E 


* 


•2a60o0/B  0_-a60o, 


-2e 


(19) 


This  expression  could  then  be  used  in  place  of  Equation  13  in  the  solution  to  the 
bond  breaking  equation,  i.e.,  the  tg  versus  aQ  relationship. 

The  fracture  of  polycrystalline  materials  involves  plastic  flow  and 
crack  propagation,  so  that  the  relationships  previously  derived  must  be  modified 
to  account  for  nonhomo geneous  rates  of  bond  rupture.  Waldorf's  material  model 
assumes  a  distribution  of  microcracks  which  grow  until  they  join  to  form  a  con¬ 
tinuous  fracture  surface.  The  lengths  and  spacing  of  the  cracks,  which  are  assumed 
to  be  nucleated  by  the  usual  mechanisms,  are  assumed  to  be  about  the  same  as  the 
grain  diameter,  so  that  each  crack  must  roughly  double  its  length  for  fracture  to 
occur.  It  is  assumed  further  that  crack  widths  are  about  equal  to  the  thickness 
of  the  grain  boundaries,  i.e.,  10  atoms. 


The  time  required  for  neighboring  cracks  to  coalesce,  and  hence  the  delay 


time  for  fracture,  t  ,  is  given  by 


t  -  - 


(20) 


where 

v  «*  the  crack  propagation  velocity  (this  is  a  function  of  the 
applied  stress) , 

b  -  one-half  the  mean  tip-to-tip  spacing  of  the  cracks. 
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In  order  to  find  the  relationship  between  crack- propagation  velocity  and 
applied  stress,  it  is  assumed  that  the  variation  of  tensile  stress  with  distance 
from  the  crack  tip  is  described  by  a  step  function.  Thus,  for  a  distance  dQ  from 
the  crack  tip,  the  stress  is  constant  and  equal  to  the  maximum  value  of  concen¬ 
trated  stress  at  the  crack  tip.  At  distances  f"om  the  crack  tip  greater  than  dQ, 
the  stress  is  constant  and  equal  to  the  applied  stress,  aQ. 


Now,  by  assuming  that  the  stress  level  in  the  region  dQ  of  stress  con¬ 
centration  is  so  large  that  fracture  will  occur  across  the  width  dQ  a  time  very 
small  compared  to  fracture  anywhere  else,  and  considering  the  grains  to  be  perfect 
homogeneous  crystals,  it  may  be  concluded  that  the  crack  tip  moves  across  the  dis¬ 
tance  d  in  the  time  t'  required  for  a  homogeneous  material  to  fracture  under  the 
concentrated  stress.  As  the  crack  propagates,  this  process  repeats,  so  that  the 
propagation  velocity  may  be  written  as 


v 


(21) 


where  t^  is  the  delay  time  for  fracture  in  homogeneous  materials  subjected  to 
constant  applied  stress. 


Expressing  dQ  and  b  in  teems  of  the  grain  diameter,  Y,  and  the  grain 
boundary  thickness,  c,  and  making  a  correction  for  plastic  flow,  the  equations  for 
propagation  velocity  and,  hence,  fracture -delay  time  may  be  obtained  in  terms  of 
the  applied  stress  and  atomic  characteristics  of  the  material.  For  a  constant 
value  of  applied  stress,  the  delay  time  versus  stress  relationship  becomec 


where 


,  a  1  (sX2^) 
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Y  «*  grain  diameter, 
c  *  grain  boundary  thickness. 
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Since  the  crack- propagation  velocity  cannot  exceed  the  velocity  of  sound 
in  the  material,  cg,  the  minimum  value  of  delay  time,  occurring  at  high  stress 
levels,  would  be  given  by 


As  a  partial  experimental  verification  of  his  theory,  Waldorf  cites  the 
agreement  between  his  predictions  and  stress-fracture  delay-time  data  obtained  by 
other  investigators  for  homogeneous  plastics  and  metals  at  relatively  low  values 
of  applied  stress  and  long  time  durations,  i.e.,  quasi-static  strengths. 

Louie  and  others  '  attempted  to  test  the  validity  of  Waldorf's  theory 
by  conducting  crack-propagation  experiments  with  plastics  and  polycrystalline 
aluminum.  By  applying  quasi-static  tensile  forces  perpendicular  to  the  line  of 
small  surface  cracks,  they  found  that  the  cracks  widened  prior  to  a  sudden  increase 
in  length.  This  circumstance  made  it  difficult  to  measure  the  crack  velocities 
accurately  and  correlate  them  with  the  applied  tensile  stresses.  Nevertheless, 
their  results  showed  no  indication  that  crack  velocity  is  an  exponential  function 
of  applied  stress,  as  predicted  by  Waldorf's  theory.  In  addition,  photographs 
showed  that  the  propagation  was  accompanied  by  considerable  plastic  flow  at  the 
crack  tips.  From  these  observations,  it  was  concluded  that  omission  of  plastic 
flow  effects  in  the  Waldorf  equation  is  a  serious  omission,  since  this  flow  can 
appreciably  alter  the  effective  stress  at  the  crack  tips  under  suitable  conditions. 
This  phenomenon  will  also  modify  the  temperature  dependence  of  the  strength  of 
the  material  on  the  &  priori  determined  parameters  of  the  theory. 

From  these  considerations,  it  would  appear  that  the  Waldorf  model  in 
its  present  form  is  applicable  only  to  very  brittle  materials.  It  should  be  noted, 
however,  that  stress  wave  propagation  in  uniaxial  strain  inhibits  plastic  flow 
due  to  the  triaxiaiity  of  the  stress  state  and  the  nigh  strain  rates.  Thus,  the 
results  of  the  crack- propagation  study  done  by  Louie  and  others  do  not  prove  con¬ 
clusively  that  the  Waldorf  model  is  unsuited  to  the  prediction  of  spall  in  ductile 
materials.  Experimental  verification  of  the  Waldorf  theory  for  the  stress  levels 
and  pulse  durations  associated  with  plate-slap  impact  spallation  has  yet  to  be 
provided. 


-33- 


SUMMARY 


Three  major  categories  of  spall  fracture  in  materials  have  been  identi¬ 
fied.  In  the  order  of  decreasing  associated  stress  levels,  they  are;  ultimate  or 
cohesive  spall,  phase  transformation  spall,  and  ductile  or  low-pressure  spall. 
Low-pressure  spall  has  been  further  subdivided  into  complete  spall,  intermediate 
spall,  and  incipient  spall.  These  low-pressure  spall  thresholds  mentioned  here 
are  also  listed  in  order  of  decreasing  associated  stress  level.  Damage  done  to 
materials  at  the  low-pressure  spall  thresholds  varies  from  separation  of  large 
parts  of  the  material  at  the  complete  spall  threshold  to  microcrack  formation  at 
the  incipient  spall  threshold. 

Several  investigators  have  found  that  other  parameters,  in  addition  to 
a  unique  value  of  tensile  stress  at  the  spall  plane,  are  necessary  to  the  formu¬ 
lation  of  an  adequate  spall  criterion.  Parameters  that  have  been  found  experi¬ 
mentally  or  offered  in  the  way  of  conjecture  as  influencing  the  spall  behavior 
of  materials  are;  tensile  pulse  duration  at  the  spall  plane,  strain  and  stress 
rates  as  the  spall  plane  is  loaded  in  tension,  spatial  stress  gradients  repre¬ 
senting  the  unloading  of  the  spall  plane  from  a  maximum  value  of  tension,  pre¬ 
compression  of  the  material  at  the  spall  plane,  ambient  temperature,  yield  strength, 
and  metallurgical  characteristics  such  as  grain  orientation. 

Spall  models  which  would  appear  to  have  some  feasibility  for  the 
prediction  of  spali  at  intermediate  stress  levels  are  neither  abundant  nor 
verified.  Four  basically  different  types  of  models  exist:  the  Cumulative  model, 
the  Energy  model,  the  Stress-Rate  model,  and  the  Bond-Breaking  model.  The  first 
three  of  these  models  are  of  a  phenomenological  or  empirical  nature,  while  the 
last  represents  an  attempt  to  describe  time-dependent  fracture  in  terms  of 
processes  taking  place  at  the  atomic  level  of  activity. 

The  cumulative  model  purports  to  predict  the  spall  due  to  an  arbitrarily 
shaped  pulse  by  combining  previously  obtained  rectangular  pulse  data  in  a  cumula¬ 
tive  manner.  Thus,  if  the  relationship  between  spall  delay  time  and  constant 
applied  tensile  stress  is  known  for  all  values  of  stress,  the  arbitrary  pulse  in 
question  can  be  approximated  by  a  series  of  rectangular  pulses  whose  time  dura¬ 
tions  can  be  compared  with  the  data.  The  ratio  of  the  duration  of  the  approximating 
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pulse  to  the  spall  delay  time  corresponding  to  that  stress  level  is  calculated 
for  each  approximating  pulse  consecutively  and  added  to  the  sum  of  those 
calculated  previously.  Spall  occurs  when  the  accumulated  sum  of  these  ratios 
is  equal  to  unity. 


The  Energy  model  conceived  at  Battelle  is  based  on  the  criterion  that 
the  mechanical  energy  supplied  to  the  spall  plane  up  to  the  time  of  fracture  is 
invariant.  Thus,  spall  would  be  predicted  to  occur  at  the  time,  t,  when  the 
following  equation  is  satisfied: 


K 


JoV."  ; 


(24) 


where 

as  m  ffs(0  “  net  tensile  stress  at  the  spall  plane, 

Vs  ”  V8(t)  *  net  particle  velocity  at  the  spall  plane, 

K  «  the  invariant  spall  energy  in  units  of  energy  per  unit  area. 


The  Stress-Rate  model  uses  the  criterion  that  the  tensile  i  trength  of 
a  material  is  dependent  on  its  tensile  unloading  characteristics,  i.e.,  stress 
rate  or  stress  gradient.  For  a  given  material,  a  locus  of  tensile  stress  versus 
corresponding  stress-rate  (or  gradient)  points,  dividing  spall  and  "no  spall", 
can  be  determined  empirically.  To  predict  spall  due  to  an  arbitrary  pulse, 
tensile  stresses  and  stress  rates  are  calculated  and  compared  with  the  critical 
values  lying  on  the  locus. 

The  Bond-Breaking  model  represents  an  attempt  to  predict  spall  from 
physical,  rather  than  phenomenological  considerations.  For  a  microscopically 
homogeneous  material,  i.e.,  a  single  crystal,  the  relationship  between  time  re¬ 
quired  for  fracture  and  applied  tensile  stress  is  found  by  integrating  an  equation 
expressing  the  rate  at  which  atomic  bonds  are  broken  as  a  function  of  stress. 

Using  this  relationship  to  determine  an  expression  for  crack- propagation  velocity 
in  terms  of  stress,  and  assuming  a  typical  size  and  distribution  of  microcracks, 
a  delay  time  versus  applied  stress  relationship  is  found  for  polycrystalline 
materials.  Since  the  delay  time  versus  applied  stress  relations  are  derived  from 
physical  considerations,  the  constants  in  the  equations  represent  physical 
characteristics  of  the  material  rather  than  the  empirical  "black  boxes"  of  the 
preceding  models.  Some  of  the  material  constants  that  are  required  are  the 
vibrational  frequency  of  the  lattice,  the  sublimation  energy,  and  grain  size. 
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APPL  iDIX  A 


SASIC  PRINCIPLES  OF  CONTINUUM  THEORY 


The  fundamental  principles  of  continuum  theory  are  restated  briefly  in 
this  Appendix.  Although  it  is  possible  to  present  these  principles  in  their 
general  tensorial  form,  this  would  only  tend  to  obscure  the  fundamental  simplic¬ 
ity  of  the  postulates.  These  are 

CONSERVATION  OF  MASS, 

CONSERVATION  OF  MOMENTUM, 

CONSERVATION  OF  MOMENT- OF-MCMENTUM, 

CONSERVATION  OF  ENERGY. 

If  thermodynamic  effects  are  included  in  a  continuum  formulation,  to  the  above 
must  be  added  the 


PRINCIPLE  OF  ENTROPY. 


In  addition,  the  basic  concepts  of 

STRESS,  STRAIN 


must  be  included. 


The  above  postulates  pertain  to  all  media  and  do  net  contain  the 
distinguishing  characteristics  of  a  particular  material.  There  must  be 
included,  therefore,  a  material  characterizing  relation  in  the  form  ot  a 

CONSTITl TIVE  EQUATION 

that  relates  the  basic  state  variables  (e.g.,  stress,  internal  energy,  and 
temperature) . 

It  is  to  tnls  basic  equation,  descriptive  of  a  material  and  its 
relation  to  the  stress  wave  propagation  problem.,  that  attention  is  directed  in 
the  following  appendices.  The  various  postulates  are  given  mathematical 
formulations  convenient  to  la  applications  at  hand. 


{BORIS' 


APPENDIX  B 


ELASTIC  WAVES 


APPENDIX  B 


ELASTIC  WAVES 


When  the  propagation  of  plane  compressional  waves  in  an  infinite, 
isotropic,  elastic  medium  Is  considered,  the  previously  stated  postulates 
yield  the  equation  of  motion: 

(X  +  2d)  ~  -  P  ~  ;  (B-l) 

dxZ  atz 


where  u  is  displacement  in  the  x  direction,  and  X,  \i  ore  materiel  moduli  (the 
Lam^  constants).  The  constitutive  equation  is  Hooke's  Law,  given  as 

°ij  ’  Xskk  6ij  +  2wij  1  (B_2) 


where  cr^  is  the  stress,  the  strain,  and  6^  the  Kronecker  delta,  and  the 
summation  convention  for  repeated  subscripts  is  utilized.  The  strain  and 
displacements  are  relate  by  the  infinitesimal  strain  expression, 
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(B-3) 


The  stated  case  of  propagation  represents  motion  under  conditions  of 
one-dimensional  strain.  Thus,  by  replacing  index  notation  with  coordinate 
subscripts,  one  has 
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The  governing  equation  of  motion  (Equation  B-l)  is  the  wave  equation,  having 
the  general  solution 


•i  =  f (x  +  c ^ t)  +  g(x  -  c^t)  ,  (B-5) 

1 

where  “[(X  +  2M>)/p]2  is  the  dilatation  velocity.  The  solution  indicates  that 
a  given  wave  shape  propagates  undistorted  with  distance  and  at  a  constant 
velocity,  c^,  as  shown  in  Figure  B-l. 
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FIGURE  B-l.  UNDISTORTED  STRESS-PULSE  PROPAGATION 


The  reflection  of  an  elastic  wave  from  a  free  boundary  results  in 
tension  waves  reflecting  as  compressive  waves,  and  incoming  compressive  waves 
reflecting  as  tension  waves.  This  latter  phenomenon  is  the  basic  mechanism 
of  spall  fracture,  both  in  the  elastic  and  inelastic  regimes. 

Although  the  character  of  elastic  waves  is  altered  by  interaction 
with  boundaries,  the  waves  do  not  interact  with  themselves.  That  is,  leftward  - 
and  rightward- propagating  elastic  stress  waves  may  encounter  and  pass  through 
one  another  without  consequent  alteration  of  their  wave  shapes  or  propagation 
directions.  During  the  duration  of  ♦'heir  encounter,  the  resulting  stress  and 
displacement  fields  will  merely  be  the  superposition  of  the  individual  fields. 
It  will  be  found  that  the  simple  reflection  and  interaction  behavior  reviewed 
here  for  elastic  waves  does  not  hold  for  elastoplastic  waves. 


Finally,  it  is  useful  to  review  propagation  of  elastic  waves  in  thin, 
straight  rods,  since  many  studies  of  dynamic  material  properties  use  this 
geometry.  The  bar  geometry,  which  is  the  simplest  case,  is  assumed  to 
represent  a  one-dimensional  stress  situation  and  has  the  governing  equation 


E 


(B-6) 


where  E  is  Young's  modulus  and  is  related  to  the  previously  given  Lam4  constant 
by  X  *  Ev/ (1  +  v)  (1  -  2v).  Only  the  c;  stress  is  nonzero,  with  all  other 

XX 

stresses  assumed  zero;  then 
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There  are  nonzero  lateral  strains,  but  since  the  inertia  associated  with  the 
resultant  motions  is  neglected,  the  remaining  strain  expressions  are  not  required. 

The  governing  equation  is  also  the  wave  equation  predicting  constancy 
of  wave  shape.  The  propagation  velocity  is  given  as  c_  **  *jE/p ,  and  is  somewhat 

D 

less  than  the  one-dimensional  strain-wave  velocity.  For  example,  in  steel, 

4  4 

Cj  *  22.4  x  10  in/sec  and  cg  *  20  x  10  in/sec.  This  decreased  velocity  for 
the  bar  case  may  be  interpreted  as  a  result  of  relaxation  of  lateral  constraints. 
For  elastic-plastic  waves,  the  plastic  wave  velocities  are  drastically  different 
between  one-dimensional  ntress  and  strain  cases,  whereas  the  velocities  given 
above  show  only  a  slight  dl.  rence  for  the  elastic  case. 


V 
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APPENDIX  C 


SHOCK  WAVES 


Ihc  preceding  discussion  of  elastic  waves  represents  wave  propagation 

4 

at  very  low  stress  regimes  (of  the  order  10  psi) *  Under  conditions  of  extremely 
high  impulsive  stress  (pressures  of  the  order  lO^psi),  shock  waves  develop  and 
propagate  in  a  material  in  a  manner  similar  to  the  fluid  dynamics  situation. 

Thus,  there  is  a  propagating  single  wave  front  across  which  the  material 
variables,  such  as  density,  stress,  velocity,  and  energy,  undergo  large  changes 
in  short  distances.  Also,  it  becomes  reasonable  to  consider  the  solid  as 
behaving  like  a  compressible  fluid,  described  by  an  equation-of-state  in  the 
form  P  *  P(v) ,  where  v  is  the  specific  volume.  Consequently,  shock  wave 
propagation,  as  in  the  case  of  elastic  waves,  has  inherent  simplifying  features 
(simplified  equation-of-state  versus  linear  elastic  behavior)  permitting  ready 
solution  to  the  stress  wave  problem,  at  least  for  the  one-dimensional  strain 
case. 

The  theory  of  shock  wave  propagation,  again  in  contrast  to  the  elastic- 
plastic  case,  has  been  reviewed  extensively,  (bee,  for  example,  References 
C-l  and  C-2.)  As  a  result,  only  one  or  two  major  aspects  of  the  theory,  useful 
for  later  comparisons  with  the  elastic-plastic  case,  will  be  reviewed  here. 

Consider  a  shock  wave  propagating,  as  shown  in  Figure  C-l. 


U  *  shock  velocity 


P  »  P 
ro’  r0 


FIGURE  C-l.  PROPAGATING  SHOCK  FRONT 


The  conservation  of  mass,  momentum,  and  energy-equations  applied  across  the 
shock  front  are,  respectively: 


P0U  *  PX(U  -  ux),  (a) 

pl  •  po  "  PoUV 


(b) 
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„  2 
P  Uu. 

»1°1  *  — 2_  +  PoU<El  •  V-  <C) 


(C-l) 


According  to  these  equations,  velocities,  densities,  pressures,  and  energies  are 
discontinuous  across  the  shock  front,  but  they  are  related  as  given  above.  There 
is  a  total  of  eight  ',>arameters  in  the  three  equations  above  (pQ,  Pp  pQ,  pp  U3 
Up  Eq,  E^).  If  it  is  assumed  that  Eq,  pQ,  and  pQ  are  known,  three  equations 
in  five  unknowns  remain.  By  eliminating  U  and  u^  between  two  of  the  equations, 
the  equation  can  be  reduced  to 

ei  -  eo  -  1/2  <  r  -  >  <pi  +  po>  •  <c-2) 

ro  rl 

This  establishes  the  states  (Ep  pp  p^)  that  can  be  reached  by  a  shock  transition 
from  an  initial  state  (Eq,  pq,  pQ).  By  elimination  of  the  particle  velocity,  Up 
from  the  mass  and  momentum  equations  above,  an  equation  for  the  shock  velocity 
becomes 

1  i 

U  »  UPl  “  P0>/(vo  *  Vj)]2  ;  (C-3) 

where  vq  =*  1/p^  and  v^  ■  1/p^  are  the  spacific  volumes. 

When  a  plane  shock  wave  arrives  at  a  free  surface,  the  shock  pressure, 

Pp  is  reduced  to  zero  by  a  rarefaction  wave  from  the  free  surface.  Both  the 
arriving  compression  and  reflected  rarefaction  waves  impart  particle  velocities 
to  the  free  surface.  The  total  free  surface  motion,  Up  is  given  by 

+  u-  ;  (c-4) 

f  c  r 


where  u£  and  uf  are  the  compression  and  rarefaction  particle  velocities, 
respectively.  Experiments  have  shown  that  U£  *  uf.  Hence,  from  Equation  C-4, 
it  is  concluded  that  the  particle  velocity,  u,  is  equal  to  one-half  the  free 
surface  velocity. 


(C-5) 


C-3 


Experimental  measurements  enable  one  to  obtain  the  shock  velocity  and 
free  surface  velocity.  These  results,  coupled  with  the  conservation  equations, 
are  sufficient  to  determine  the  equation-of-state  for  the  material.  As  an 
example,  consider  the  results  obtained  from  an  experiment  on  Armco  iron.(^”^ 
Measurements  indicated  that 

Pq  »  7.87  g/cc, 

U  =  6.06  x  10^  cm/sec, 
uf  **  2.786  x  105  cm/sec, 

Pq  *»  1  x  10  dynes/cm  .  (C-6) 


From  Equation  C-5, 


u  =  1.393  x  1(T  cm/ sec  . 


By  substituting  in  Equation  C-l  (b)  ana  solving  for  p^,  ore  obtains 

9  2 

p^  «*  665  x  10  dynes /cm  =  665  kbar  . 


(C-7) 


(C-8) 


To  obtain  the  specific  volume  v^,  the  above  data  is  substituted  in  Equation 
C-3  and  solved  by  v^,  giving 

Vj  «*  0.09787  gram/cm^  .  (C-9) 

Consequently,  v^/vq  =  0.77,  indicating  that  the  iron  was  compressed  to  77%  of  its 
original  volume. 

The  results  of  a  series  of  experiments  at  differing  shock  pressures 
constitute  a  Hugoniot  equation-of-state  curve  which  relates  pressure  to  specific 
volume.  A  typical  result  for  copper  is  shown  in  Figure  C-2. 

Because  of  permanent  deformations  and  heat  gen«'ration  that  occurs 
during  the  shock  transition,  the  Hugoniot  is  neither  an  adiabatic  r.or  an 
isothermal  curve.  It  is,  instead,  a  dynamic  compressibility  curve  representing 
an  irreversible  process;  that  is,  -  Sq  >  0.  However,  more  detailed  thermo¬ 
dynamic  considerations  show  that  dS  is  small,  and  that  it  is  a  reasonable 


C-4 


approximation  to  consider  the  process  as  isentropic  (i.e.,  constant  entropy), 
This  approximation,  which  is  useful  at  shock  pressure,  will  also  be  used  in 
the  elastic-plastic  wave  process. 
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FIGURE  C-2.  HUGONIOT  CURVE  FOR  COPPER^"3) 
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APPENDIX  D 


ELASTIC- PLASTIC  WAVE  PROPAGATION 

The  past  sections  have  reviewed  briefly  the  basic  principles  of  con¬ 
tinuum  theory  and  some  of  the  major  features  of  elastic  and  shock  wave  propaga¬ 
tion.  The  cases  of  wave  propagation  considered  represented  extremes  in  terms  of 
stress  levels,  with  elastic  waves  corresponding  to  wave  propagation  at  low  stress 
levels  and  shock  waves  corresponding  to  extremely  high  stress  levels.  It  was 
noted  that  at  either  extreme,  various  assumptions  tending  to  simplify  the  propa¬ 
gation  problem  became  reasonable.  The  present  case  of  elastic-plastic  waves 
represents  the  intermediate  case  of  stress  level. 

At  intermediate  stress  levels,  the  inherent  ability  of  solid  materials 
to  support  significant  shear  stresses  makes  the  assumption  of  ideal  fluid 
behavior  used  in  shock  analysis  a  less  accurate  one.  Thus,  neglecting  material 
rigidity,  which  represents  a  realistic  assumption  at  extremely  high  stresses, 
becomes  less  plausible  at  intermediate  stress  levels  and  must  be  replaced  by 
a  more  accurate  model  of  the  material. 

In  this  section,  the  theory  of  plasticity  will  be  applied  to  the 
derivation  of  the  governing  equations  for  stress  wave  propagation  under  con¬ 
ditions  of  one-dimensional  strain.  A  number  of  complicating  subsidiary  effects, 
known  to  exist  in  some  materials,  will  be  neglected  in  this  preliminary  review 
of  theory,  but  will  be  considered  in  later  sections.  Attention  will,  instead, 
be  focused  on  the  major  features  of  wave  propagation  and  the  governing  equations. 


The  Longitudinal  Stress-Strain  Rela tionshii 


In  this  section,  the  essential  features  of  stress-strain  curves 
resulting  from  qvcei-static  >  uniaxial**  stress  test  conditions  will  be  reviewed, 
as  well  as  the  idealized  forms  used  for  mathematical  analysis.  The  construction 
of  the  elastic-plastic  stress-strain  relations  corresponding  to  conditions  of 
uniaxial  strain  will  then  be  given. 


D-2 


Stress-Strain  Relations  Under  Uniaxial  Stress 


Representative  results  of  stress-strain  curves  obtained  for  various 
materials  under  quasi-static,  isothermal,  uniaxial  tensile- stress  conditions 
are  shown  in  Figure  D-l. 
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FIGURE  D-l.  TYPICAL  STRESS-STRAIN  CURVES  FOR  VARIOUS  MATERIALS 


The  curves  are  not  meant  to  represent  any  specific  material,  but,  instead,  to 
illustrate  a  number  of  properties.  Thus,  material  (a)  is  a  linear  elastic 
material  until  the  proportional  limit,  Op,  is  reached.  Beyond  ap,  the  material 
deforms  plastically.  The  increasing  stress  with  strain,  however,  indicates  that 
work-hardening  is  occurring.  Material  (b)  possesses  a  definite  yield  point, 
ov.  In  some  cases,  a  lower  yield  point  exists,  as  indicated  by  the  dotted  line. 
Material  (c)  has  no  well-defined,  proportional  limit.  Of  course,  such  material 
behavior  represents  the  macroscopic  response  of  the  material  and  is  a  consequence 
of  the  microscopic  deformation  processes  reviewed  earlier  in  chis  report. 

The  behavior  described  above  is  a  result  of  a  continually  increasing 
load.  If  the  specimen  is  unloaded,  reverse  loaded,  or  load  cycled  during  a  test, 
responses  as  depicted  in  Figure  D-2  are  possible. 


FIGURE  D-2 


TYPICAL  LOAD  CYCLING  AND  REVERSE  LOADING  CURVE 


D-3 


Thus,  the  material  has  been  loaded  beyond  its  yield  point,  and  then  unloaded  (a). 
Upon  reloading,  a  small  hysteresis  loop  is  formed.  Both  unloading  and  loading  are 
elastic  until,  in  the  latter  case,  the  load  reaches  the  original  unloading  value, 
(b) ,  at  which  yielding  again  occurs.  If,  at  a  later  stage,  unloading  and  then 
reverse  loading  occurs,  reverse  yielding  takes  place  approximately  at  (c) ,  which 
is  the  magnitude  of  the  maximum  stress  attained  less  twice  the  yield  stress,  a y. 
However,  reverse  yielding  may  occur  at  a  somewhat  lesser  value,  (d) ,  in  which 
case  the  material  is  said  to  exhibit  the  Bauschinger  effect. 

A  detailed  mathematical  representation  of  the  actual  stress-strain 
curve  is  generally  too  complicated  for  use  in  analysis.  Hence,  simplifications 
are  often  adopted  which,  although  approximating  the  actual  stress-strain  rela¬ 
tions,  greatly  reduce  the  complexity  of  the  mathematical  representation.  Figure 
D-3  shows  a  number  of  idealized  curves. 
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FIGURE  D-3.  MATHEMATICAL  IDEALIZATIONS  OF  STRESS-STRAIN  CURVES 

In  Figure  D-3,  (a)  is  a  bilinear  representation  of  an  elastic-plastic 
material  with  work  hardening;  (b)  is  an  elastic,  perfectly-plastic  material. 

The  lower  curve  is  an  idealization  of  several  aspects  of  unloading  and  reverse 
loading.  Thus,  the  small  hysteresis  loop  which  occurs  during  load  cycling  has 
been  neglected,  (c) .  Reverse  yielding  is  shown  at  (d) ,  and  a  representation  of 
the  Bauschinger  effect  is  indicated  by  (e) . 
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The  stress-strain  relations,  discussed  in  the  foregoing,  have  been 
based  on  engineering  stress  and  strain  defined  as 


L  -  L 


G  =  A  ,  6 


(D-l) 


However,  these  definitions  are  often  replaced  by  true  stress  and  true  strain 
definitions,  given  as 

-  P  _  ,  L  . 

<?»—,£  -  In  ”  > 


(D-2) 


where  A  represents  the  instantaneous  cross-sectional  area,  and  L  is  the  current 
length.  The  relationships  between  the  true  stresses  and  strains  and  the  engineer¬ 
ing  counterparts*  are  given  by 


YTe  »  G  *  ln^  +  e>  • 


(D-3) 


Stress-Strain  Relations  for  Uniaxial  Strain 

The  governing  equations  of  elastic-plastic  wave  propagation  for  con¬ 
ditions  of  one-dimensional  strain  require  the  appropriate  stress-strain 
relations.  These  will  be  derived  in  this  section  using  plasticity  theory  and  the 
quasi-static,  stress-strain  relations  from  uniaxial  stress  conditions.  The  basic 
concepts  of  plasticity  theory  will  be  reviewed  briefly  before  developing  the 
stress-strain  relations  using  the  Fowles  procedure. 

Plasticity  Theory.  When  a  material  is  loaded,  permanent  deformation 
occurs  as  stresses  exceed  certain  limits  characteristic  of  the  material,  including 
its  past  loading  history.  Unloading  allows  a  degree  of  elastic  recovery.  In 
plasticity  theory,  a  basic  assumption  is  made  that  there  exists  a  scalar 

p 

function,  f(<?ijt  eij>  H)»  caHed  a  yield  function,  which  depends  on  the  state 
of  stress  and  strain  and  the  history  of  loading.  This  function  characterizes 


*  In  the  work  that  follows,  the  engineering  definition  is  used. 
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the  yielding  of  the  material  under  any  possible  combination  of  stresses.  Here, 

p 

is  the  stress  state,  e^j  is  the  plastic  strain,  and  k  is  a  parameter 
representing  the  previous  loading  history.  The  equation,  f  ■  0,  represents  a 
closed  surface  in  the  stress  space.  No  change  in  plastic  deformation  occurs 
when  f  <  0.  Plastic  deformation  only  occurs  for  f  ■  0.  No  meaning  is  associated 
with  f  >  0,  provided  the  material  properties  are  independent  of  strain  rate. 

Yield  functions  of  complete  generality,  containing  the  parameters 
indicated  above  and  successfully  describing  all  aspects  of  plastic  deformation 
of  materials,  have  not  been  developed.  However,  many  materials  may  be 
adequately  described  by  simplified  forms  of  the  yield  function.  Thus,  if  the 
plastic  deformation  characteristics  of  a  material  are  isotropic  and  work¬ 
hardening  effects  are  neglected,  it  is  reasonable  to  assume  that  plastic  yielding 
can  depend  only  on  the  magnitudes  of  the  three  principal  applied  stresses,  o^, 
a2*  aV  and  not  on  t*ietr  directions.  Hence,  a  yield  function  of  the  form^D'2^ 


f(J1?J,,,J3)  -  0 


may  be  assumed,  where  J^,  is  the  invariant  in  each  of  the  first  three  stress 
tensors  given  by 

+  o2  +  Oy  J2  -  -0^2  +  a2<73  + 


J3  "  Ws  * 


(D-4) 


(D-5) 


However,  experimental  observations  have  indicated  that  plastic  yielding  is  quite 
insensitive  to  hydrostatic  stress.  By  introducing  the  deviatoric  principal 
stresses,  ct^,  Oy  0'3>  defined  as 


a'  *  a  -  —  (a  +  a  +al 

1  1  3  '12  y  ’ 

*  n 

az  x  °2  '  3  +  a2  +  ’  °3  *  °3  3  ^al  *  C2  +  °3)  ’ 


(D-6) 


and  noting  that  *^i  *  ai  +  a2  °3  *  0>  the  general  form  of  the  plastic  yield 
criteria  may  be  reduced  to 


fCJj.  J3)  -  0  . 


(D-7) 
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Many  specific  forms  of  yield  criteria  have  been  suggested  over  the 
years.  However,  most  of  these  have  attributed  some  influence  to  hydrostatic 
stress  and,  thu<  ,  have  not  been  in  accord  with  experimental  observations.  The 
two  simplest  criteria  not  having  this  fault  were  those  suggested  by  Tresca 
(1864)  ’"d  von  Mi  .'es  (1913).  The  Tresca  criterion  states  that  yielding  occurs 
when  the  maximum  shear  stress  reaches  a  certain  value;  thus, 

o^  -  o^  *•  constant  (D-8) 

where  >  cr^  >  ° y  The  above  may  be  expressed  in  terms  of  J^,  Jy  but  the 
results  are  somewhat  cumbersome. (D“3) 

The  von  Mises  criterion  does  not  involve  the  function  J^,  but  states 
that  yielding  occurs  when  J reaches  a  critical  value,  given  as 

+  a'£  +  ol^  “  constant;  (D-9) 


where  the  deviatoric  stress  tensor,  ct1^  ,  *8  defined  in  terms  of  by 


o’  m  o  -  06 
ij  ij  ij  * 


(D-10) 


Both  models  describe  elastic,  perfeetly-plastic  behavior.  For  most  metals,  it 
is  found  that  the  von  Mises  criterion  fits  experimental  data  more  closely,  although 
the  Tresca  condition  is  simpler  to  uce  in  theoretical  applications .  This  diffi¬ 
culty  is  sometimes  resolved  by  using  the  Tresca  condition  with  an  empirical 
adjustment  factor  intended  to  minimize  the  differences  between  Lie  two  results. 
Fortunately,  in  the  one-dimeneional  situations  to  be  considered,  both  criteria 
reduce  to  the  same  result. 

The  two  major  mathematical  formulations  of  plaeticity  are  identified 
as  plastic  deformation  theory  (also  known  as  total  strain  theory)  and  plastic 
flow  theory  (also  known  as  incremental  theory).  Although  the  deformation  theory 
finds  application  in  some  problems  involving  continuous  loading,  it  yields 
anomalous  predictions  under  conditions  of  cyclic  or  reverse  loading.  The  p-astie 
flow  theory,  which  rc,'  tes  increments  of  strain  to  increments  of  stress  at  any 
given  state  of  strets  and  strain,  has  generally  found  the  widest  acceptance  and 
application  and  is  utilized  in  the  present  review, 
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In  the  elastic  region,  the  relation  between  stress  and  strain  is 
given  by  Hooke’s  Law.  In  incremental  tensor  form,  this  is 


do’ 


de 


ij 


ii 


2>a 


2V 


6 .  do 

ij 


(D-1X) 


where  do 


ij 


is  the  deviatoric  stress  increment  defined  as 


dcf^  *  do  -  do  ,  (D-12) 

and  do  *  do^  /3  is  the  hydrostatic  component.  Compressive  stresses  and  strains 
are  taken  as  positive. 

When  the  material  is  strained  beyond  the  yield  point,  the  total 
strain  increment,  de^at  a  given  state  of  strain  is  comprised  of  an  elastic 
and  a  plastic  portion  given  as 

de^  »  ds^j  +  de^j  (D-13) 


Hence,  ac  a  given  stress-strain  state, 
be  given  as 


the  increment  of  plastic  work,  dW  ,  will 

P 


d«p  -  0lj  d.' 


The  von  Mises  yield  criterion  has  the  form 

’ij 


V<«  )  -  (f)l  (o'  o i(J  , 


(D-14) 


(0-15) 


where  Y  is  the  yield  stress  in  simple  tension.  Since  Y  may  change  as  a  result 

of  work  hardening  (i.e.,  plastic  work,  Wp) ,  this  is  accounted  for  by  the 

functional  form  Y  (W  ) ,  and  is  illustrated  in  figure 

P 
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FIGURE  D-4.  CHANGE  OF  YIELD  POINT  AS  A  RESULT  OF  PLASTIC  WORK 

Yj  *  initial  yield  point 

Y  «*  increasing  value  of  yield  point 
as  a  result  of  plastic  work. 


Finally,  it  was  previously  pointed  out  that  plastic  deformation  was 
observed  to  be  insensitive  to  hydrostatic  pressure.  This  observation  is  stated 
mathematically  by  requiring  no  permanent  volumetric  change  due  to  plastic 
deformation;  it  is  given  as 

de^  «  (D-16) 

The  fundamental  relations  of  elastic-plastic  theory,  mentioned  above,  will  now 
be  specialized  to  one-dimensional  stress  and  one-dimensional  strain  situations. 


One-Dimensional  Stress.  The  cne-dimeneional  stress  situation  corre¬ 
sponding  to  uniaxial  tension  and  compression  test  conditions  is  approximately 
achieved  during  longitudinal  wave  propagation  in  rods.  This  phenomenon  will  be 
developed  briefly  here.  To  distinguish  between  later  one-dimensional  strain 
results,  the  direction  of  stress  will  be  taken  here  as  the  s  direction.  Then, 
the  incremental  elastic  stress-strain  relations,  Equation  D-ll,  become 


de 


e 

s 


3E 


da 


(D-17) 


The  von  Mises  yield  criterion  (Equation  D-15)  reduces  to  c  -  Y(W  ) „ 

*  '  P 

The  expression  for  the  increment  of  plastic  work,  dWn,  reduces  to 

d.W  *  Y  dep  ; 
p  s 


(D-18) 
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where  dep  *  de  -  de6.  If  there  is  no  work  hardening,  Y  ■  constant,  and  dY  *»  0 
s  s  s 

for  the  added  strain  increment.  Consequently,  de  *  0.  For  the  more  general  case 

s 

of  work  hardening  (i.e.,  dY  f  0) ,  dee  ■  dY/E  so  that  Equation  D-16  becomes 

8 

dW  -  Y[de  -  da  /E]  -  Y[d$  -  dY/E].  (D-19) 

p  8  S  8 

The  overall  relationship  of  total  work,  plastic  work,  and  increments  of  total 
work  and  plastic  work  is  illustrated  in  Figure  D-5.  A  diagram  that  helps  clarify 
the  remarks  on  the  effects  of  work  hardening  on  plastic  strain  and  that  led  to 
Equation  D-19  is  given  in  Figure  D-6. 


FIGURE  D-5.  WORK  DONE  DURING  PLASTIC  DEFORMATION 

(a)  W  *  plastic  work,  W  *  elastic  work,  and  W  +  W  =  total  work  at  stress  level  o, 

p  e  p  e 

(b)  For  added  strain  de  at  stress  a,  dW  ■  total  added  work,  dW  *  ode  ■ 

irrecoverable  plastic  work.  P  P 


FIGURE  D-6.  EFFECTS  OF  WORK  HARDENING  ON  THE  RELATIONSHIP 
BETWEEN  TOTAL  ELASTIC  AND  PLASTIC  STRAIN 

For  work  hardening,  de^  *  deP^,  dee  *  aY/E. 

0  p 

If  work  hardening  is  absent,  dY  a  0,  de  =  0,  de^  =  de  ^  =  dt . 
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One-Dlmensicmal  Strain.  The  x  direction  will  be  designated  as  the 
strain  direction.  The  elastic  relations.  Equation  D-ll,  then  become 

d6x  *  E  (dJx  "  2V  dV’  d6y  "  d6z  “  E  (1  "  V)day  "  Vdax  5  (D'20) 

where  the  equality  do  “  da  ,  due  to  symmetry,  has  been  utilized  in  the  above. 

y  z 

The  plastic  incompressibility  relation,  Equation  D-16,  becomes 

dep  +  2dcp  -  0  .  (D-21) 

x  y 

The  von  Mises  criterion  (Equation  D-15)  reduces  to 

a  -a  -  Y(W  )  .  (D-22) 

x  y  p 

It  should  be  recalled  that  the  von  Mises  and  Tresca  yield  criteria  are  identical 
under  the  conditions  of  uniaxial  stress  and  uniaxial  strain,  as  above. 

Finally,  the  general  relation  for  plastic  work  must  be  specialized  to 
the  present  case.  Thus,  Equation  D-14  reduces  to 

dW  -  a  dep  +  2o  dep  .  (D~23) 

p  x  x  y  y 

However,  dep  «  -dep/2  froin  Equation  D-21,  so  that  the  above  becomes 

dW  «  (a  -  a  )deF  .  (D-24) 

p  x  y  x 

Using  the  yield  condition,  Equation  D-22  thus  becomes 

dW  -  Y(W  )d€p  ,  (D-25) 

p  p  x  ’ 

which  is  the  same  as  the  case  for  one- dimensional  stress. 

In  order  to  express  the  incremental  plastic  work  in  terms  of  total 
strain  increment  and  work-hardening  characteristics  in  a  manner  analogous  to 
Equation  D-19  for  one-dimensional  stress,  the  following  relation  is  utilized: 


(D-26) 


D-ll 


Recall  that,  for  the  one-dimensional  stress  situation,  the  case  of  no  work 

6  6 

hardening  yielded  de  ■  0;  whereas  for  work  hardening.de  *  dY/E  was  obtained. 

8  8 

In  the  present  case,  from  Equation  D-ll, 


.  e  x  A  (l-2v)  _ 

dSx  ■  d0  ! 


(D-27) 


so  that 


da’ 

deP  ■  dc - - 

x  X  2(1 


;i-2v) 


(D-28) 


Expressions  may  be  developed  for  da^  ,  da  that  make  it  possible  to  put  the  above 
in  the  desired  form.  Thus  it  is  observed  that 


d€ii  '  deii  +  deiim  dsii 


(D-29) 


by  virtue  of  Equation  D-16.  However,  *  dex’  8^nce  de  *  ^ez  “ 

e  p  Y 

since  de  *  de  -  de  ,  Equation  D-2S  may  be  written  as 

X  X 


de  =dee  .  !lii  + 
d  ii  dii  2p  E 


(D-30) 


Noting  that  da^  ,  do^  (*  da^)  are  increments  of  the  deviatoric  principal  stresses, 
and  recalling  that  in  an  earlier  observation  =  ai  +  a£  +  a^  =  0,  one  finds 
that  da^  *  0.  Consequently, 


deli  ■  dtx 


CD— 31 > 


It  is  further  noted  that 


da'  «=  da  -  da  «  da  -  \  (da  +  2da  )  =  \  (da  -  da  )  .  (D-32) 

xx  x  3  x  y  3  x  y 


From  the  yield  criterion.  Equation  D-22,  da  -  da  «*  dY  so  that  da'  =  —  dY, 

x  y  x  j 

Substituting  these  results,  and  substituting  Equation  R-31  in  Equation  D-28,  one 


obtains 


.  p  2  ,  .  dY  , 
d€x  =  3  ld€x  ‘  *  ]  • 


(D-33) 
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This  permits  expression  of  Equation  D-25  in  the  desired  form: 


Y  ‘de*  -  S  »  • 


(D-34) 


The  next  objective  is  to  relate  explicitly  the  stress  increment,  da  , 

6  P 

to  the  total  strain  increment,  de  .  Now,  de  »  de  +  der  .  The  conditions  of 

X  X  X  X 

p  6 

uniaxial  strain  require  that  dey  ■  dez  *  0,  so  that  de^  *  -  dey  .  Using  this 
and  Equation  D-21,  deP  *  -  2deP  “  2  dc6  ,  so  that  the  total  strain  increment, 

x  y  y 

dex»  is  given  by 


de  “  dg6  +  2deS 
x  x  y 


(D-35) 


^  6 

Using  Hooke's  Law,  Equation  D-ll,  for  de  ,  de  ,  one  obtains 

x  y 


de. 


(1-2*)- 

E 


(dax  +  2day)  . 


(D-36) 


Now,  the  yield  criterion  may  be  put  in  the  form 


da  +  2 da  -  3da  -  2dY. 
x  y  x 


(D-37) 


Hence,  Equation  D-36  becomes 


de  «  ILLiv)  -  |  dY) 

x  E  x  3 


(D-3S) 


Prom  the  relationships  between  the  elastic  moduli,  the  bulk  modulus,  K,  is  given 
by  K  -  E/3(l-2v)  so  that 

da  -  K  de  +  \  dY  ,  (D-39) 

X  X  J 


which  is  the  desired  :  ;remental  stress-strain  relation.  The  mean  stress,  da, 

2 

is  given  by  da  »  Kde  •*  Kde  .  Hence,  Equation  D-39  becomes  da  -  da  *  -r  dY. 

X  X 

This  integrates  to 


a  -  a 
x 


-•  y. 
3 


I 

I 

I 


(D-40) 
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Thus,  beyond  the  material  yield  limit,  the  difference  between  tne  one-dimensional 

strains  and  hydrostatic  stress-strain  curves  is  only  a  function  of  the  material 

yield  stress  as  a  function  of  strain.  If  the  material  does  not  strain  harden, 

2 

dY  *  0,  Y  ■  Yo  *  constant,  and  a  -  o  *  —  Yo. 

x  3 

Some  of  the  features  of  the  resulting  stress-strain  relations  for 
one- dimensional  strain  are  reviewed  in  Figure  D-7. 


FIGURE  D-7.  STRESS-STRAIN  RELATION  FOR 
ONE-DIMENSIONAL  STRAIN 

(a)  Material  hydrostat;  da^/de^  *  K  *  slope. 

(b)  Perfectly-plastic  material;  parallel  to 
hydrostat  but  removed  a  distance  Y  . 

(c)  Work-hardening  material;  dcr^  =  Kde^  +  "f 


r~* 


\ 


Li  t 

j 
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and  they  are  nearly  the  same  for  most  materials.  However,  beyond  the  yield 
point,  the  slopes  differ  radically.  Thus,  for  a  non-work-hardening  material, 
the  slope  of  the  one-dimensional  strain  case  is  da^/de^  «  K,  whereas, 
dcr  /de  ■  0.  This  factor  will  be  of  significance  when  wave-propagation 

8  S 

velocities  are  considered. 

The  preceding  development  has  served  to  delineate  the  major  features 
of  the  stress-strain  curve,  such  as  the  slopes  of  the  curve  (before  and  after 
yielding)  and  the  relation  to  the  hydrostat.  It  is  evident  that  the  parameters 
necessary  for  the  construction  of  the  curve  are  the  elastic  constants  and  the 
yield  stress  obtained  from  uniaxial  stress  experiments.  In  fact,  if  the  stress- 
strain  curve  for  uniaxial  strain  is  to  be  approximated  by  straight  line  segments 
with  an  abrupt  break  at  the  yield  stress,  the  preceding  is  sufficient.  However, 
if  a  more  detailed  picture  of  the  transition  region  between  elastic  and  plastic 
behavior  is  desired,  further  considerations  are  necessary. 

It  is  possible  to  obtain  the  desired  information  from  uniaxial  stress 
data,  provided  the  comparison  is  made  at  corresponding  values  of  strain  for 
which  the  plastic  work  is  the  same.  Equating  increments  of  work,  dW^,  for  the 
two  cases,  as  given  by  Equations  D-19  and  D-34,  gives 

f  -  f 1  ’  Y[de„  -  f 1 '  •  (D'41> 

This  simplifies  to 

de*  -  1  des  ■  i  •  (D-42) 


If  any  variation  in  K  with  strain  is  neglected,  the  above  may  be  integrated  to 
give 


e 

x 


o 

e 

x 


Slzfl 

6K 


> 


where  the  "o"  superscript  refers  to  initial  yield  values.  However, 


1  +v  v°  c 

Y  ’  es 


,o 


(D-43) 


o 


E 


(D-44) 
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so  that  Equation  D-43  reduces  to 

ex  *  2  es  '  6K  *  <D-45) 

Thus,  the  procedure  is  to  calculate  e  from  uniaxial  stress  data, 
using  Equation  D-45,  and  to  obtain  cr^,  using  Equation  D-40  and  the  fact  that 
<y  *  Ke^.  Figure  D-9  illustrates  the  procedure. 


FIGURE  D-9.  RELATION  OF  UNIAXIAL  STRAIN 
AND  STRESS  CURVES 


A  slightly  modified  procedure,  described  by  Barker,  et  al,(®“4)  calculates  the 
strain  offsets  relative  to  the  elastic  slopes  in  order  to  establish  the  com¬ 
parison. 

At  this  stage,  a  constitutive  relation  between  stress  and  strain  for 
uniaxial  strain  conditions  has  been  derived.  The  uniaxial  stress,  isothermal, 
quasi-static  stress-strain  curve  has  been  utilized  in  the  development.  Varia¬ 
tions  in  the  elastic  constants  due  to  finite  strain  effects  have  been  neglected. 


The  applicability  of  the  resulting  quasi-static  relation  to  the  dyn'unic  situation 
of  elastic-plastic  wave  propagation  is,  of  course,  open  to  question.  Ultimately, 
the  degree  of  comparison  between  experimental  results  and  theoretical  predictions 
will  establish  the  applicability.  These  will  be  discussed  in  a  later  section  in 
an  attempt  to  arrive  at  an  assessment  of  such  items  as  strain  rate,  finite  strain, 
thermodynamic  and  Bauschinger  effects.  For  the  immediate  purpose  of  reviewing 
wave  propagation,  unobscured  by  the  complicating  effects  mentioned  above,  the 
quasi-static,  isothermal  relation,  as  given  by  Equation  D-39,  will  be  used. 


The  Basic  Governing  Equations 


The  basic  postulates  of  continuum  mechanics  will  be  applied  to  derive 
the  governing  equations  of  elastic-plastic  wave  propagation.  The  constitutive 
relations  developed  in  the  preceding  section  will  be  utilized  and  the  method 
cf  characteristics  will  be  applied  to  solve  the  equations.  The  development 
parallels  that  of  Morland . 

Wave  propagation  in  the  x  direction  will  be  assumed,  and  a  Lagrangian 
coordinate  system  (x,t)  will  be  employed.  Under  these  conditions,  conser¬ 
vation  of  mass  requires  that 

P (1  -  ex)  =  ;  (D-46) 


where  p,  and  P  are  initial  and  final  densities,  respectively;  e  is  the 
i  X 

infinitesimal  compressive  strain  defined  as  =  -  du/dx;  and  u(x,t)  is  the 
particle  displacement.  Conservation  of  momentum  gives 


Sv 

at 


(D-47) 


where  v(x,t)  is  the  particle  velocity  defined  as  v  =  3u/St.  The  compressive 

stress  is  a  .  Because  uniaxial  strain  conditions  will  be  assumed,  there  will 

also  exist  stresses  a  and  o  ,  where  o  =  o 

y  z  y  z 
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Assuming  the  existence  of  a  constitutive  equation,  a  *  a  (e  ) ,  as 

X  XX 

discussed  in  the  previous  section,  the  above  may  be  written  as 


dt 


2  9ex 
(T(e  )~ 
X  ox 


(D-48) 


where 


o  i  da 

C  (e  )  «  —  ~ — 

x  p,  de 

i  x 


From  the  definitions  of  strain  and  velocity, 


5v 

5x 


(D~49) 


(D-50) 


Using  Equation  D-50,  Equation  D-48  may  also  be  written  as 


+  C2(e  ) 
x 


(D-51) 


The  moment- of -momentum  postulate  of  continuum  mechanics  does  not  enter 
here  because  of  the  uniform  normal  stresses  and  zero  shear  stresses.  Similarly, 
the  energy  equation  does  not  directly  influence  the  development.  Neglecting 
thermodynamic  effects,  it  would  merely  represent  a  first  integral  of  the 
momentum  equation. 

The  method  of  characteristics  will  be  used  to  solve  the  system 
governing  equations,  given  by  Equations  D-48  and  D-50.  The  characteristics  in 
the  x  -  t  plane  for  these  equations  are  given  by 

dx/dt  -  +  C(ex),  or  ~  -  +  1/C(ex>  .  (D-52) 


The  invariant  quantities  along  these  characteristics  are 
,  e 

“  +  C(e  ),  v  x  f  x  C(e  )de  »  constant  »  R,  (D-53) 

and 

.  e 

dt  “  '  v  "  $  X  C^Sx^dex  *  con8Cant  “  R2*  (D-54) 

where  the  constants  (Rj,  R0)  are,  in  general,  different  along  different 
characteristics . 
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Further 
fellows:  Suppose 
material,  whare  v 
before  arrival  of 
region  of  the  x  - 


information  regarding  the  above  constants  may  be  developed  as 
the  stress  pulse  is  propagating  into  a  uniform  region  of  the 
and  are  constant.  This  would  correspond  to  the  situation 
the  pulse,  where  v  *  «  0.  It  is  represented  by  the  shaded 

t  plane  in  Figure  0-10. 


FIGURE  D-10 .  THE  x  -  t  PLANE,  INCLUDING  A  UNIFORM 

REGION  OF  v  -  e  -  0 
x 


Since  each  negative  characteristic  passes  through  the  uniform  region,  it  may  be 
concluded  that  the  invariant,  R£,  is  the  same  for  all  characteristics  and  may  be 
replaced  by  R.  A  positive  characteristic  is  intersected  at  each  point  by  a 
negative  characteristic,  so  that  along  each  positive  characteristic 


e 

v+r  c<ex)dex 


€ 

R. ,  v  -  T  x  C(e  }de  *>  R. 

1  «  X  X 
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Results  of  adding  and  subtracting  the  above  ate 


Rx  +  R 


J  11  cep*. 


v  “ 


2 


2 


(D-56) 
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e 

T.n  other  words,  v  and  J  xC(ex)dex  aru  separately  constant  along  each  positive 
characteristic.  Further,  since  C(e  )>0,  J  xC(e  )de  is  an  increasing  function 
of  e  ,  it  is  concluded  that  c  must  remain  constant  along  a  positive  character- 

X  X 

istic  if  the  integral  is  to  remain  constant.  It  follows  directly  that  the 
dependent  variables,  C(e  ),  a  (e  ),  are  also  constant.  A  further  result, 

X  XX 

useful  in  a  later,  limiting  cases  is  obtained  by  integrating  Equatica  D-47. 
Thus, 
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v 
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Propagation  of  Elastic-Plastic  Waves  -  General  Description 

The  propagation  of  an  elastic-plastic  wave  considered  here  is  that 
initiated  by  a  general,  prescribed  stress  applied  on  the  surface  of  a  half¬ 
space.  The  essential  feacures  of  the  initial  stages  of  propagation  will  be  the 
existence  of  an  elastic  precursor  or  forerunner,  a  propagating  plastic  wave 
front  lagging  behind  the  elastic  front,  an  elastic  unloading  front  that  is 
overtaking  the  plastic  front,  and  finally,  a  possible  plastic  unloading  front. 

The  formation  of  plastic  shock  fronts  will  be  considered,  as  well  as  specialized 
loading  pulses.  (However,  treatment  of  complicated  wave  interactions  that  occur 
when  unloading  waves  overtake  the  loading  fronts  will  be  treated  in  the  next 
section.)  The  equations  of  motion  and  constitutive  relations  developed  in  the 
previous  section  will  be  utilized.  Thus,  complications  due  fco  strain-rate  or 
Bauschinger  effects  will  be  avoided  at  this  stage,  as  well  as  other  subsidiary 
effects  briefly  discussed  in  the  previous  sections. 

The  prescribed  stress  pulse  is  shown  in  Figure  D-ll(a).  Figure  D-ll(b) 
shows  the  associated  stress-strain  curve.  Each  stress  amplitude  propagates  with 
the  local  slope  of  the  stress-strain  curve.  The  associated  Lagrangian  diagram 
is  shown  in  Figure  D-12. 
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The  portion  OA  of  the  pulse  propagates  with  the  elastic  velocity, 

Cq.  These  are  represented  on  the  Lagrangian  diagram  by  the  rays  emanating 
between  0  and  t^,  the  time  duration  of  the  elastic  portion  of  the  pulse.  The 
AB  portion  of  the  pulse  is  represented  by  rays  from  t^  -  t^  of  th»  Lagrangian 
diagram.  The  various  stress  levels  of  thi3  portion  propagate  with  different 
velocities,  as  indicated  by  the  variation  of  slope  of  the  stress-strain  diagram 
in  the  AB  region.  Thus,  because  of  the  upward  curvature  of  AB,  higher  values 
of  stress  propagate  at  higher  velocities.  This  is  indicated  by  the  clustering 
together  of  the  rays  from  tA  -  tg.  A  region  of  constant  stress,  BC,  then 
exists  from  tg  to  t^,. 

Unloading  occurs  along  CE.  The  initial  stage,  CD,  is  elastic  and  is 
indicated  in  the  Lagrangian  diagram  by  rays  from  tc  -  The  final  stage  of 

unloading,  DE,  is  anelastic  due  co  reverse  yielding  and  has  various  velocities 
associated  with  the  various  stress  levels.  The  rays  from  t^  -  t£ 
represent  this  region.  The  divergence  of  this  family  of  rays  is  a  consequence 
of  unloading  from  higher  stress  with  associated  higher  velocities  to  lower 
stress  with  associated  lower  velocity  levels. 

It  is  of  interest  to  consider  the  variations  of  pulse  shape  that  arise 
from  variations  of  wave  velocity  indicated  in  the  foregoing.  There  would,  of 
course,  be  no  such  variations  if  the  loading  were  within  the  elastic  region, 
since  the  stress-strain  curve  would  have  constant  slope.  All  rays  of  the 
Lagrangian  diagram  would  then  be  parallel  and  would  indicate  undistorted 
propagation  of  the  stress  pulse.  However,  variations  do  occur  for  the 
anelas  ic  case,  as  shown  in  Figure  D-13,  where  the  pulse  shapes  at  two  time 
instants  are  shown.  The  major  features  of  figure  D-13  are:  (a)  the  elastic 
precursor  or  forerunner,  OA;  (b)  the  plastic  wave,  AB,  whose  front  is 
increasingly  steep  due  to  the  greater  propagation  velocity  of  the  higher  stress 
levels;  (c)  the  region  of  constant  stress,  BC,  that  is  diminishing  because  of 
the  overtaking  elastic  unloading  wave,  CD;  (d)  the  plastic  tail,  DE,  that  is 
becoming  stretched  out  because  of  the  slower  propagation  velocity  of  the  lower 
stresses.  Because  the  effects  described  are  essentially  related  to  the  stress- 
amplitude  dependence  of  the  velocity,  the  subsequent  wave  distortion  is  referred 
to  as  amplitude  dispersion. 
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FIGURE  D-13.  VARIATIONS  IN  STRESS  PULSE  AT  TWO  INSTANTS  OF  TIME 


This  description  of  the  progressive  variation  in  the  wave  3hape  has 
been  terminated  before  two  major  events  have  occurred.  These  are  the  formation 
of  a  plastic  shock  front  at  AB  and  the  overtaking  of  the  loading  front  (either 
before  or  after  shock  formation)  by  the  unloading  front,  CD.  Consideration  will 
be  given  to  shock  formation.  The  stages  leading  to  this  situation  are  shown 
in  Figure  D-14.  The  steepening  of  the  front  ie  shown  for  t  =  tQ,  t^, 


FIGURE  D-14.  FORMATION  OF  THE  PLASTIC  SHOCK  FRONT 


t 


At  t  =  ty  a  physically  unacceptable  situation  has  arisen  as  the  more  rapidly 
traveling  stress  components  have  "overlapped”  the  slcver,  lower  stress  level 
portion  of  the  wave.  This  situation  is  indicated  on  the  Lagrsngian  diagram 
by  the  crossing  of  characteristics. 

It  is  when  the  wavefront,  or  a  portion  of  the  wavefront,  becomes 
vertical  that  the  continuous  plastic  wave  front  breaks  down  and  a  discontinuous 
shock  front  is  formed  traveling  at  a  shock  velocity,  U.  Across  the  shock  front, 
the  variables  of  stress,  density,  velocity,  and  internal  energy  vary  discon- 
tinuously  as  shown  in  Figure  D-15. 
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FIGURE  D-15.  SHOCK  FRONT  ACROSS  WHICH  MATERIAL 
VARIABLES  CHANGE  DISCONTINUOUSLY 


The  equations  of  conservation  of  mass,  momentum,  and  energy  still 
pertain  in  the  region  separating  the  two  states,  o and  If  q^  *  vq  -  U  and 

*  v^  -  U  define  the  material  velocities  relative  to  the  shock,  the  mass  and 
momentum  equations  become 


Plql  3  Poqo  *  CT1  "  G0  3  poqo2  -  plql2 


(D-58) 


Those  may  be  solved  to  give 


pi<°l  -  ao) 


Po<°l  ' 


•  po>'  1  pI(pl  -  po>  ' 


(D-59) 
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Conservation  of  energy  is  expressed  by 


"l’l*!  '  Wo  '  "Vo  -  al’l  +  2  Poqo3  '  2  Ol’l3  -  <D-60> 


where  Eq,  are  the  internal  energies.  When  the  values  of  q^,  qQ  given  in 
Equation  D-59  are  substituted  in  the  energy  expression,  one  obtains 
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This  is  the  well-known  Rankine-Hugoniot  relation  presented  earlier  in  the  review 
of  hydrodynamic  shock  theory.  By  substituting  for  p  in  terms  of  E,  the  results 


for  q  ,  q,  become 


-d-€o) 


Gl-°o 
pi(ereo) 


-i  I 

2 


(1-ep 


o  -a 
1  o 

pi(ere0) 


(D-62) 


Finally,  by  defining  the  particle  jump  velocity  across  the  shock  as  v 


v,  -  v 
1  o 
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To  complete  the  development,  it  would  be  useful  to  compare  the  pre¬ 
dicted  velocity  jump  via  the  shock  mechanism  to  that  predicted  by  the  continuous 
solution.  The  details  of  this,  which  require  assumption  of  the  stress-strain 
curve  form,  will  not  be  presented  here.  However,  Mori  and  (^5)  has  investigated 
this  and  found  that  the  difference  between  the  two  solutions  is  only  about  0.37.. 
The  magnitude  of  the  shock  velocity,  U,  is  given  by 
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Of  particular  interest  in  Equation  D-64  is  the  form  of  the  square  root  term, 
which  can  be  interpreted  in  terms  of  the  chord  connecting  the  two  points  on  the 
stress  strain  curve  connected  by  the  shock  transition.  This  is  shown  in 
Figure  D-16. 


The  location  of  the  shock  front  in  space  and  time  must  be  considered. 
That  is,  the  location,  x,  where  the  shock  front  forms  must  be  established,  as 
well  as  the  subsequent  location  of  the  front,  xg(t).  The  details  of  this,  which 
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shock  Chat  occurs  for  a  smooth  loading  now  occurs  instantaneously  under  step 
loading.  This  is  illustrated  in  (b)  by  the  shock  chord  connecting  and 
and  in  (c)  by  only  two  rays  from  the  origin.  The  resulting  wave  shape  is  shown 
in  (d) .  Fortunately,  the  simplified  case  discussed  is  the  one  of  greatest 
interest  since  it  represents  the  nature  of  the  stress  pulse  induced  in  plate- 
impact  experiments. 

Before  concluding  this  section  on  shock  formation,  it  is  of  interest 
to  note  other  forms  the  elastic-plastic  wave  front  may  take  due  to  variations 
in  the  stress-strain  curve  of  the  material.  These  are  shown  in  Figure  D-18, 
along  with  wave  fronts  characteristic  of  such  materials. 


FIGURE  D-18.  VARIATIONS  IN  WAVEFORM  DUE  TO  MATERIAL  CHARACTERISTICS (D-6) 

S- 

(a)  Material  without  a  well-defined  yield  point. 

(b)  Material  without  a  linear  elastic  portion. 


Elastic-Plastic  Waves  of  Unloading  and  Reflection 

The  interaction  between  loading  and  unloading  portions  of  the  pulse 
occurs  when  the  unloading  wave  overtakes  the  loading  wave  front.  At  this  point, 
it  is  encountering  a  change  in  acoustic  impedance  and  reflects  from  and  transmits 
through  this  discontinuity  accordingly.  The  general  features  of  the  reflection 
of  elastic-plastic  waves  from  a  free  surface  will  also  be  described. 


D-28 


The  approach  taken  by  Morland  is  to  approximate  the  wave  profile  by 
a  series  of  discontinuous  steps,  or  stress  increments,  and  to  consider  the 
interaction  of  these  increments.  A  stress  increment,  6a,  propagating  at  a 
velocity,  c,  will  cause  a  change  of  particle  velocity,  6v,  and  strain,  fie. 

These  quantities,  which  result  from  a  discontinuous  pulse,  may  be  calculated 
from  Equation  D-57,  which  was  developed  for  the  continuous  case  by  considering 
the  limiting  case  of  the  pulse  shown  in  Figure  D-19. 


^0 
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FIGURE  D-19.  PROPAGATING  STRESS  PULSE 
OF  LENGTH  fix 


Thus,  in  Equation  D-57,  it  may  be  written  that 


da  a  -  a 
_ x  1  o 

dx  fix 
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Noting  that  fix  =  cfit,  it  can  be  found  that 
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The  jump  in  strain  is  found  from  the  definition  of  c,  i.e.,  c 


fie  ■  fia/P^c 


fia/pjfie, 


Thus, 

(D-67) 


Now  consider  the  situation  of  an  elastic  unloading  stress  increment 
overtaking  a  plastic  loading  increment,  as  shown  in  Figure  D-20.  The  dotted  lines 
represent  the  point  of  interaction;  or,  P  are  particles  on  either  side  of  this  line. 
After  interaction,  a  new  stress,  o^,  exists  at  the  interaction  line.  This  may  be 
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FIGURE  D-20.  SITUATION  BEFORE,  (a)  AND  AFTER  (b)  INTERACTION 
BE1VEEN  ELASTIC- PLASTIC  STRESS  INCREMENTS <D~ 5 > 

determined  by  applying  velocity  continuity  conditions  subsequent  to  the 

interaction.  Thus,  it  may  be  written  for  v^,  that 

a2  -  °i  c  -  o2  a3  -  aQ  a  -  a 

a  1  picl  Pico  pico  P  1  picl 
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The  final  velocity  of  O'  is  determined  by  the  initial  velocity,  v^,  and  the 
velocity  increments  imparted  by  the  plastic  loading  wave  Ccr2~cri^  >  the  elastic 
unloading  wave  (aQ“a^) »  and  the  wave  resulting  from  the  interaction  (cr^-a^) • 
For  only  the  latter  of  these  three  waves  affects  its  initial  velocity,"  v^. 
Writing  v^  »  gives 


c  -  c  2.C 

o  -  o  -  — - -  A’  ,  a_  -  a.  =  A - - - A'  ; 

3  o  c  +  c.  3  1  c+c. 
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where  A, A'  are  the  loading  and  unloading  stress  increments,  respectively,  given 
by 


A  “  °2  "  ai  ’  A'  =  CT2  "  ao  ’ 


(D-70) 


The  resulting  strains,  6^,6^,  may  be  calculated  by  using  the  incre¬ 
mental  expression,  Equation  D-67.  Thus, 
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From  Equations  D-69  and  D-70,  this  gives 
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This  indicates  that  a  discontinuity  in  strain  has  been  introduced  at  the 
interaction  line,  and  that  it  subsequently  propagates  at  the  local  particle 
velocity. 


Before  the  wave  interaction  process  can  be  completed,  the  interaction 
of  the  reflected  elastic  loading  wave,  creat.  y  the  foregoing  process ,  and 
other  overtaking  elastic  unloading  waves  rout.  .e  considered.  The  situation  is 
shown  in  Figure  D-21. 
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FIGURE  D-21.  SITUATION  BEFORE  (a)  AN.,  AFTER  (b)  INTERACTION 
BETWEEN  ELASTIC  STRESS  INCREMENTS  <0-5) 


The  particle  velocities,  v^,  v^,  are  given  by 
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Putting  va  *  vg  gives 


C3  "  °o  “  C1  ~  C2  “  A  ’  (D-75) 

a  -  a  n  a.  -  j  =  A’  ; 
i  j  /  o 

where  A, A'  are  the  loading  and  unloading  increments.  Calculation  of  in  the 

manner  used  for  the  plastic  interaction,  shovB  that  *  e^,  so  that  no  strain 
discontinuity  arises. 

Thus,  if  a  series  of  elastic  unloading  increments  propagating  at 

velocity  cq  overtakes  a  plastic  front  propagating  at  c^,  the  sequence  of  events 

is  as  follows;  (1)  the  first  unloading  increment  interacts  with  the  plastic 

front;  (2)  a  reduced  plastic  front  is  formed  of  stress  and  propagates  ahead 

at  c^;  (3)  an  elastic  loading  increment  of  stress  is  reflected  at  a  velocity, 

c  ;  (4)  the  reflected  increment  interacts  with  an  oncoming  elastic  unloading 
o 

increment;  (5)  the  two  elastic  increments  pass  through  each  other,  stepped  to 
new  levels;  (6)  the  second  elastic  unloading  increment,  modified  by  the  inter¬ 
action,  overtakes  the  plastic  front  aud  interacts  with  it,  starting  the  sequence 
of  events  over  again.  Further  consideration  will  be  given  to  this  process 
in  regard  to  specific  aspects  of  the  plate  impact. 

Consideration  is  now  extended  to  the  reflection  of  waves  from  a  free 
surface.  For  simplicity,  an  elastic-plastic  wave  having  a  two-front  structure, 
as  shown  in  Figure  D-22  (a),  will  be  considered  to  impinge  on  the  free  surface. 

By  simple  superposition  methods,  it  is  established  that  the  elastic  compression 
front  reflects  from  the  surface  as  a  tension  front  and  induces  a  change  of 
velocity.  The  reflected  tension,  propagating  to  the  right  and  annihilating 
the  compressive  elastic  precursor,  effectively  creates  an  unloading  wave 
propagating  to  the  right,  as  shown  in  (b)  of  the  figure.  The  elastic  unloading 
wave  then  interacts  with  the  oncoming  plastic  front,  reduces  the  magnitude  of 
the  plastic  wave,  and  creatt *  a  reflected  elao'zic  loading  wave  propagating 
back  to  the  free  surface  ahead  of  the  reduced  plastic  front  (c) . 
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FIGURE  D-22 .  VARIOUS  STAGES  OF  WAVE  REFLECTION  AND 
INTERACTION  WITH  A  FREE  SURFACE 


It  is,  of  course,  evident  that  the  above  interaction  process  could 
continue  an  indefinite  number  of  times  before  the  plastic  front  reaches  the 
surface.  In  actuality,  the  amplitude  of  the  reverberating  elastic  wave  is 
rapidly  reduced  to  an  insignificant  amount  compared  to  the  plastic  stress 
amplitude,  and  only  the  first  reverberation  is  of  significance.  Thus,  the 
first  reflection  of  the  elastic  front  from  the  plastic  wave  will,  when  it 
reaches  the  free  surface,  induce  an  additional  step  in  velocity.  The  plastic 
wave  front  will  then  reach  the  surface,  stepping  the  velocity  yet  again.  The 
velocity  and  displacement  profiles  of  the  free  surface  motion  would  have  the 
general  appearance  shown  in  Figur  :  D-23.  Further  consideration  of  the  actual 
velocity  changes  induced  will  be  given  in  a  subsequent  section  in  which  the 
entire  wave  propagation  process  due  to  plate  impact  is  considered. 
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FIGURE  D-23.  VELOCITY  AND  DISPLACEMENT  PROFILES 
OF  FREE  SURFACE  MOTION 


Plate  Impact  and  Determination  of  Material  Properties 

In  the  last  three  sections,  various  aspects  of  elastic-plastic  wave 
initiation,  propagation,  interaction,  and  reflection  have  been  considered. 

In  the  present  section,  the  wave  propagation  problem  resulting  from  the  planar 
impact  of  two  flat  plates  will  be  reviewed.  For  such  a  problem,  conditions 
of  uniaxial  strain  pertain  in  the  interior  regions  of  the  plates  until  the 
first  rarefactions  from  the  edges  reach  the  center.  All  aspects  of  initiation, 
propagation,  etc.,  discussed  previously,  will  be  present.  Although  the 
complexities  of  t.ie  many  wave  interactions  will  preclude  extensive  hand 
analysis,  the  early  stages  of  the  process  will  be  reviewed.  A  relatively  simple 
elastic,  perfectly-plaetic  material  model  will  bp  used  so  that  the  basic 
features  of  the  stress  wave  system  may  be  easily  studied.  The  method  of 
deducing  the  stress-strain  relation  of  the  material  from  the  measurement  of 
rear  surface  motion  will  be  reviewed. 

An  example,  taken  from  Penning,  et  al,^”^)  wiH  be  used  to  illustrate 
various  aspects  of  the  resultant  stress  wave  system.  The  planar  impact  of  a 
driver  plate  upon  a  target  plate,  where  both  plates  are  of  the  same  material, 
will  be  observed.  An  elastic,  per fectly-plastic  material,  having  the  charac¬ 
teristics  mentioned  and  illustrated  on  pages  D-l  through  D-17  and  Figure  D-8, 
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will  be  used  to  represent  the  material.  The  Lagrangian  wave  diagram  will  be 
employed  to  depict  the  wave  process. 

Example.  Consider  the  case  of  a  driver  plate  impacting  the  target 
plate  at  a  velocity  of  0.010  cm/psec,  inducing  a  compressive  stress  of 
8.1  kbar  at  the  plate  interface.  This  compression  will  be  achieved  in 
two  steps.  There  will  be  a  compression  to  6.33  kbar  propagating  into 
the  target  a*  the  elastic  velocity  of  0.623  cm/p-sec,  followed  by  a  further 
compression  to  8.1  kbar  propagating  into  the  target  at  0.541  cm/nsec. 

The  former  wave  is  the  elastic  precursor,  and  the  second  wave  represents 
the  plastic  front.  Exactly  the  same  thing  is  occurring  in  the  driver 
plate.  This  stage  of  the  process  is  indicated  by  OA  and  OP  on  the  wave 
diagram  of  Figure  D-24. 

The  compressive  waves  first  reach  a  free  surface  at  P  in  the  driver 
plate.  The  elastic  front  reflects  as  a  decompression,  or  unloading 
wave,  releasing  the  material  from  6.33  kbar  to  zero  stress.  This 
front  interacts  with  the  oncoming  compressive  front  of  magnitude 
1.77  kbar  at  point  Q.  The  analysis  of  the  interaction  at  Q  follows 
the  methods  outlined  on  pages  D-27  through  D-33.  The  resulting  stress 
wave  system  is  composed  of  a  recompression  or  loading  wave  of 
amplitude  1.96  kbar,  propagating  on  to  the  free  surface  at  the 
elastic  velocity,  and  an  unloading  wave  from  3.1  kbar  to  1.96  kbar, 
propagating  (at  the  elastic  velocity)  back  through  the  driver  plate. 

At  the  interaction  line  between  these  two  waves,  a  density  discon¬ 
tinuity  or  contact  surface  has  been  created.  The  recompression  wave 
reflects  from  the  driver- free  surface  as  an  unloading  wave  and 
propagates  back  through  the  driver  plate  leaving  unstressed,  unstrained 
material  m  its  wake  until  it  reaches  the  contact  surface.  Beyond  the 
contact  surface,  the  material  is  left  unstressed,  but  with  a  net 
residual  strain.  Exactly  the  same  interaction  occurs  in  the  target  plate 
at  A,  only  at  a  later  time.  Lwo  stages  of  this  interaction  process  are 
shown  in  Figure  D-25.  One  occurs  at  t  ■  t^,  subsequent  to  the  reflection 
from  the  driver-free  surface,  but  prior  to  reflection  from  the  target 
surface.  The  other  is  at  t  =  t^,  subsequent  to  the  reflection  from  the 


FIGURE  D-24.  WAVE  DIAGRAM  FOR  AN  8.1  KBAR  COMPRESSION  PULSE(D~7) 
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target  surface,  but  prior  to  further  interaction.  These  times  are 
indicated  on  the  Lagrangian  diagram,  Figure  D-24,  as  dashed  lines 
at  t  *  tj.t^. 

At  this  stage,  the  interaction  between  the  unloading  waves  from 
the  driver  and  target-free  surface  begins.  The  first  interaction  is 
between  the  8.1  kbar  to  1.96  kbar  unloading  waves,  which  occurs  at 
point  B  of  the  Lagrangian  diagram.  This  produces  two  oppositely-moving 
tension  waves  of  magnitude  -4.28  kbar  propagating  at  the  elastic 
velocity.  These  interact  almost  immediately  with  the  oncoming  1.96 
kbar  to  0  unloading  fronts  (points  C  and  D  on  the  diagram) .  Each  of 
these  interactions  (G  and  D)  produces  two  oppositely-directed  tension 
waves  (this  is  illustrated  for  the  interaction  at  C  by  the  symbol  t_*.) . 
Now,  the  tension  waves  produced  at  C  are  actually  dual  fronts;  the  first 
step  takes  the  pressure  from  -4.28  kbar  to  -4.52  kbar  (and  propagates 
at  the  elastic  velocity  of  0.623  cm/psec),  while  the  second  step  is 
from  -4.52  kbar  to  -5.96  kbar  and  propagates  at  the  plastic  velocity 
0.541  cm/psec.  The  identical  situation  is  created  at  D.  Lack  of  space 
on  the  Lagrangian  diagram  prevents  showing  all  of  these  details.  Thus, 
of  the  two  dual  systems  from  C,  only  the  one  propagating  to  the  left 
shows  both  wave  fronts  (those  propagating  to  the  right  from  D) .  The 
resultant  stress  wave  system  is  shown  at  an  intermediate  time,  t  *  t^, 
in  Figure  D-25.  This  is  a  time  instant  between  the  identical  inter¬ 
actions  at  C  and  D  described  above. 

Now,  two  of  the  dual  systems  created  at  C  and  D  undergo  a  head-on 
interaction  at  E  of  the  diagram.  This  produces  a  region  of  -7.60  kbar 
stress  bounded  by  waves  propagating  at  the  plastic  velocity  in  the 
opposite  directions.  The  resulting  stress  situation  is  shown  at  t  *  t^ 
in  Figure  D-25. 

The  next  step  of  the  process  occurs  when  the  dua)  wave  system, 
emanating  from  D  toward  the  target- free  surface,  encounters  the  contact 
surface  formed  from  the  first  wave  reflection  off  the  free  surface 
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(identical  contact  surfaces  were  formed  in  the  driver  and  target 
plates  from  the  reflections  and  interactions  at  P  and  A).  The  first 
step  of  the  dual  wave  passes  through  at  unchanged  magnitude  and  speed. 

When  the  second  step  encounters  the  surface,  a  very  weak  elastic 
tension  wave  is  reflected,  moving  to  the  left  at  0.623  cm/psec  and 
dropping  the  pressure  about  0,15  kbar.  This  wave  is  neglected  on 
the  diagram.  The  second  step  also  refracts,  taking  on  the  elastic 
velocity  and  an  amplitude  of  -6.06  kbar.  This  situation  is  shown  as 
t  "  tj  in  Figure  D-25. 

The  next  stage  of  the  process  concerns  the  reverberation  of  the 
-4.52  kbar  tensile  wave  between  the  target-free  surface  and  the  -6.06 
kbar  tensile  wave  front.  Thus,  the  former  tensile  wave  reflects  from 
the  free  surface  as  a  -4.52  kbar  to  0  compressive  wave  at  unchanged 
velocity.  It  interacts  with  the  second  tensile  wave  step  producing  a 
compression  wave,  from  -6.06  to  -1.51  kbar  moving  toward  the  driver 
plate,  and  a  tension  wave,  from  0  to  -1.51  kbar  moving  toward  the  free 
surface  where  it  reflects  as  a  compression  wave  from  -1.51  to  0  kbar. 

This  interaction  occurs  m  the  vicinity  of  point  F  on  the  Lagrangian 
diagram.  Because  of  the  fine  detail  of  the  wave  structure  at  this 
stage,  it  is  not  possible  to  show  all  of  the  rays  of  the  diagram.  Thus, 
a  single  ray,  marked  G  in  the  diagram,  has  been  used  to  depict  the  0 
to  -1.51  and  -1.51  to  -6.06  fronts.  The  situation  is  shown  as  t  =  tg 
in  Figure  D-25. 

Meanwhile,  the  main  tension  wave  (marked  as  "main  T-wave"  on  the 
Lagrangian  diagram)  encounters  the  previously  mentioned  contact  surface. 

In  reality,  this  wave  system,  which  resulted  from  the  previously  described 
wave  interaction  at  E,  is  a  dual  wave  structure,  although  this  detail  is 
omitted  on  the  Lagrangian  diagram  (thus,  only  a  single  ray  is  shown 
propagating  to  the  right  from  E) .  At  this  stage,  the  interactions  that 
occur  with  the  contact  surfaces  and  with  the  compression  fronts  are  too 
finely  spaced  on  the  Lagrangian  diagram,  and  too  numerous  to  describe  in 
detail.  The  contact  surface  interaction  is  followed  by  wave  interaction 
at  point  H  of  the  diagram. 
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After  the  interaction  at  H,  waves  of  tensile  unloading,  marked  I 
in  the  diagram,  propagate  toward  the  driver,  while  waves  of  tensile 
loading  propagate  toward  the  target- free  surface  and  reflect.  The 
situation,  after  the  above  interactions,  is  shown  at  time  t  =  t^  in 
Figure  D-25.  Sonu  of  the  fine  details  of  the  wave  structure  in  this 
diagram  are  only  suggested  in  approximate  terms. 

The  wave  system,  for  the  example  just  considered,  was  somewhat 
simplified  by  virtue  of  the  magnitude  of  the  induced  stress  at  collision.  Thus, 
the  impact  velocity  was  such  that  the  peak  stress  was  8.1  kbar.  At  this 
level,  the  unloading  that  occurs  when  the  first  waves  reach  the  free  surfaces 
is  purely  elastic  (points  A  and  P  of  Figure  D-24).  If  the  impact  velocity 
is  high  enough,  reverse  yielding  occurs  at  unloading,  creating  an  additional 
wave  front.  This  factor  is  illustrated  in  Figure  D-26  for  a  14  kbar  impact 
stress.  At  point  P  of  that  figure,  the  plastic  front  reflects  as  an  elastic 
and  a  plastic  front,  as  indicated  by  the  dual  waves  emanating  from  that  point. 

A  similar  situation  holds  at  point  A  of  the  diagram.  The  subsequent  wave 
interactions  are  considerably  more  numerous  than  those  described  in  the 
lower  stress  level  impact  example.  If  the  impact  stress  is  still  higher,  the 
consequent  wave  interactions  are  even  mere  complicated  due,  in  general,  to  the 
reverse  yielding  effects  which  occur  during  unloading.  This  is  indicated  in 
Figure  D-27  for  a  19  kbar  impact  stress. 

The  preceding  discussion  has  been  devoted  to  the  stress  history  in  the 
interior  of  the  material  as  a  result  of  plate  impact.  The  results  have  been 
based  on  an  a  priori  knowledge  of  the  constitutive  equation  of  the  material  and 
the  conditions  of  impact.  The  inverse  problem,  in  which  the  material  properties 
are  to  be  found  frou  the  plate-impact  test,  is  of  equal  interest.  Thus,  assuming 
that  both  impact  conditions  and  experimental  data  on  the  ensuing  motion  of  the 
target-free  surface  are  known,  determination  of  the  material  constitutive 
equation  is  desired.  The  procedure  for  this  is  outlined  in  the  following. 
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FIGURE  D-26.  WAVE  DIAGRAM  FOR  A  14-XBAR  COMPRESSION  PULSE^0-7^ 


FIGURE  D-27.  WAVE  DIAGRAM  FOR  A  19-KBAR  COMPRESSION  PULSE^0"'' 
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A  typical  free-surface  motion  record  from  a  plate -impact  experiment 
is  shown  in  Figure  D-28.  From  the  discussion  of  the  last  few  pages,  it  is 


Time  After  Impact,  ftsec 

FIGURE  D-28.  FREE  SURFACE  DISTANCE-TIME  PLOT  FOR  , 
TYPICAL  PLATE-IMPACT  EXPERIMENT  (D'8 


realized  that  the  observed  free-surface  motion  is  a  result  of  a  rather  compli¬ 
cated  system  of  wave  interactions.  A  procedure  for  obtaining  the  constitutive 
equation  for  a  material  from  free-surface  motion  measurements  would  be  to 
assume  a  constitutive  equation  and  then,  using  this  relation  and  the  impact 
conditions,  to  solve  the  wave- propagation  problem.  Such  a  solution,  using  a 
characteristics  computer  code,  would  account  for  all  of  the  complex  wave 
interactions  and  would  compute  a  predicted  free-surface  motion.  The  results 
would  be  compared  with  the  experimental  data.  If  the  comparison  were  poor,  the 
constitutive  relation  would  be  revised  and  the  analysis  repeated  in  an 
iterative  manner  until  agreement  was  reached.  However,  if  the  initial  assumed 
constitutive  relation  greatly  deviates  from  the  actual  one,  a  number  of  itera¬ 
tions  might  be  required  beiore  good  match  between  theory  and  experiment  is 
obtained.  A  procedure,  which  is  described  below,  yields  a  very  good  approxi¬ 
mation  to  the  constitutive  relation  by  neglecting  many  of  the  complex  wave 
interactions. 
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The  experimentally  obtained  displacement-time  record  is  approximated 
by  a  series  of  straight  line  segments.  When  a  step-pressure  pulse  arrives  at  a 
free  surface,  it  causes  a  step  change  in  surface  velocity  which  translates  as  a 
sloped  straight  line  in  the  displacement-time  record  of  the  motion.  If  the 
reflection  and  subsequent  interaction  of  the  pressure  step  with  che  remainder 
of  the  oncoming  wave  is  neglected,  then  approximation  of  the  free- surface  data 
by  straight  line  segments  is  equivalent  to  representing  the  wave  shape  by  a 
coi'r  "  series  of  discrete  pressure  steps,  or  stress  jumps.  The  situation 
is  i 1 in  Figure  D-29,  which  is  intended  to  depict  a  stress  wave  at  an 
arbiv  •  ,  t  *  t^,  prior  to  impinging  on  the  free  surface.  Each  stress 

increaeu  -<s  considered  to  propagate  with  its  own  characteristic  velocity. 


FIGURE  D-29.  APPROXIMATION  OF  A  STRESS  PULSE  BY  STRESS  JUMPS 


The  procedure  is  to  apply  the  Rankine-Hugoniot  conditions  across  the 
stress  jumps.  Thus,  the  relations  of  interest  are  the  momentum  and  mass 
equations  which  have  the  form 
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(D-76) 


where  Ao  ie  the  stress  jump,  U  the  shock  velocity,  Au  the  change  in  particle 
X  8  P 

velocity,  and  p^,p£  densities  ahead  of  and  behind  the  shock,  respectively. 

The  parameters  determined  from  the  experiment  are  Ug,  Au^,  and  p^.  Thus, 
velocity  of  propagation  of  each  stress  step  is  found  merely  by  dividing  the  plate 
thickness  by  the  time  required  for  propagation  from  the  impact  surface  to  the 
free  surface.  The  change  in  particle  velocity  is  found  by  taking  one  half  of 
the  difference  between  successive  velocity  seeps  of  the  free  surface.  The 
density  ahead  of  the  shock,  is  considered  as  known.  Thus,  A *nd  may 
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be  found,  which  establishes  points  on  the  stres8-3train  curve. 

An  outline  of  a  typical  computation  would  be  as  follows. 
Given :  Target- plate  thickness  =  T 
Initial  target  density  =  pQ 
Displacement-time  data  given  by  Figure  D-30. 


Time 


FIGURE  D-30.  HYPOTHETICAL  DISPLACEMENT-TIME 
HISTORY  OF  REAR  SURFACE  MOTION 


Calculations : 

(1)  Velocity  of  first  shock, 

Change  of  particle  velocity  , 

Change  of  stress, 

Change  of  density, 

Resultant  strain, 

(2)  Velocity  of  second  shock, 


U, 


o 


Au, 


Aax  -  ax  -  oo  =  ot  -  PCVV 

Pj  Ui 

P0=U1-  Aul’ 
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Aa2  *  °2  “  al  "  P1U2Au2  ’ 
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In  this  way,  successive  pairs  of  points,  o^,  e^,  a^,  62*  cr^,  e^,  etc.,  are 
established  on  the  stress-strain  curve  and  connected  by  straight  line  segments. 

Having  obtained  an  approximate  stress-strain  curve  in  this  way, 
the  results  are  then  used  in  the  computer  program  which  takes  into  account 
the  complex  wave  interactions  to  predict  a  free-surface  motion.  Although  a 
number  of  approximations  are  used  in  generating  the  initial  stress-strain 
curve,  it  is  frequently  found  that  the  subsequent  exact  analysis  yields 
predicted  free-surface  motion  that  agrees,  within  experimental  error,  with 
measured  data  and  that  no  further  iterations  are  required. 

Nevertheless,  it  must  be  realized  that  the  foregoing  approximate 
procedure  neglects  density  changes  insofar  as  wave  interactions  with  these 
changes  are  concerned.  Also,  unloading  by  rarefaction  waves  has  been  neglected, 
and  the  particle  velocity  has  been  approximated  by  one-half  of  the  free-surface 
motion.  This  last,  which  holds  for  elastic  waves,  is  not  strictly  true  for  the 
inelastic  case. 
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MATERIAL  PROPERTIES  EFFECTS 

The  propagation  of  elastic-plastic  waves  under  conditions  of  one¬ 
dimensional  strain  has  been  considered  in  the  past  several  sections.  The 
equations  governing  the  propagation  (developed  on  pp  D-l-D-17)  were  derived  on  the 
basis  of  certain  kinematical  assumptions  on  the  deformation  of  the  material  and  on 
a  rather  simple  constitutive  equation.  For  example,  temperature  and  strain-rate 
effects  were  assumed  to  be  negligible.  From  the  standpoint  of  understanding  the 
essential  features  of  elastic-plastic  wave  propagation,  it  was  certainly  advisable 
to  make  such  assumptions  so  that  subtle  material  effects  would  not  obscure  the 
basic  features  of  the  propagation. 

However,  since  the  prediction  of  spallation  relies  on  determination  of 
the  stress  field  in  the  interior  of  the  material,  it  is  advisable  to  reexamine  the 
basic  governing  equations  of  the  wave  propagation  since  these  are  utilized  to  deter¬ 
mine  the  stress  field.  It  is  important,  for  example,  to  establish  whether  the 
assumed  constitutive  relation  accurately  reflects  the  real  behavior  of  the  material 
and  to  assess  the  degree  of  confidence  that  may  be  placed  in  experimentally  deter¬ 
mined  material  relations.  It  is  also  important  to  review  the  current  assessment 
of  the  value  of  strain-rate  and  temperature  effects  in  describing  the  dynamic 
response  of  materials. 


Review  of  Basic  Equations 

In  order  to  provide  a  basis  for  the  reexamination  of  the  governing 
equations,  the  equations  on  pages  D-17  -  D-20  are  restated. 

Conservation  of  mass: 

P(1  -  ex)  -  pt 

Conservation  of  momentum: 


Kinematics  (strain); 


Constitutive  relation: 


ox  “  CTx  (ex)  . 
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On  the  basis  of  elastic* plastic  theory,  the  constitutive  relation  for  a  material 
under  conditions  of  one-dimensional  strain  was  developed  on  pages  D-l  -  D-17. 

Also  illustrated  there  were  the  forms  taken  by  an  elastic,  perfectly-plastic  and 
an  elastic,  work-hardening  material.  How  to  establish  the  quasi-static  consti¬ 
tutive  relation  directly  from  the  results  of  quasi-static  uniaxial  stress  tests 
was  shown  as  well. 

The  order  of  the  reexamination  shall  now  be  as  follows: 

(a)  Thermodynamic  effects 

(b)  Work-hardening  and  Bauschinger  effects 

(c)  Strain-rate  effects 

(d)  Finite  strain  effects. 

Thermodynamic  Considerations 

Before  discussing  possible  temperature  effects  on  elastic-plastic  wave 
propagation,  it  is  desirable  to  review  some  general  thermodynamic  considerations 
of  material  stress-strain  relations.  In  a  propagating  stress-wave  system,  the 
first  physical  action  is  generally  that  which  compresses  the  material,  causing 
a  temperature  rise  of  some  magnitude.  This  action  is  followed  by  decompression 
caused  by  rarefaction  waves.  The  decompression  cools  the  material.  If  the 
process  occurs  rapidly  enough,  so  that  heat  conduction  processes  do  not  have 
sufficient  time  to  act,  and  if  initial  and  final  temperatures  are  the  same 
indicating  that  a  reversible  action  has  occurred,  the  elastic-plastic  shock-wave 
system  would  be  thermodynamically  described  as  isentropic-adiabatic.  The 
resulting  material  stress-strain  curve  would  be  an  isentropic  relation. 

However,  it  might  be  recalled  that  the  initial  development  of  the  elastic- 
plastic  stress-strain  relations  for  conditions  of  one-dimensional  strain 
(pages  D-l  -  D-17)  were  based  on  utilizing  the  results  of  quasi-static,  isothermal, 
uniaxial  stress  tests.  The  resulting  stress-strain  curve  for  uniaxial  strain  con¬ 
ditions  was,  thus,  an  isothermal  relation.  Since  the  conditions  of  stress-wave 
propagation  are  more  likely  to  be  adiabatic  conditions,  a  governing  stress-strain 
relation  based  on  isothermal  conditions  would  ’  >  in  error  thermodynamically. 

The  above  error  may  be  eliminated  by  utilizing  adiabatic  data  to  develop 
the  uniaxial  strain  relation.  Thus,  for  materials  adequately  described  by  a  simple 
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elastic,  perfectly^plastic  material  model,  the  resulting  stress-strain  relation 

is  determined  by  two  straight  lines;  one  the  elastic  slope  and  the  other  paralleling 

the  hydrostat.  The  slopes  of  these  lines  are  directly  related  to  the  elastic 

constants.  Hence,  by  establishing  these  comtants  under  adialatic  conditions  using 

ultrasonic  pulse-echo  techniques,  an  adiabatic,  uniaxial  strain  relation  may  be 

(E-l) 

constructed.  This  method  has  been  utilized  by  Barker,  et  al,  and  others. 

The  plate-slap  experimental  technique  offers  the  possibility  of  estab¬ 
lishing  the  material  constitutive  relation  directly,  in  contrast  to  conversion  of 
uniaxial  stress  test  data.  To  a  certain  extent,  consequently,  the  preceding  remarks 
on  isothermal  versus  adiabatic  conditions  must  be  reinterpreted.  The  question  now 
is  not  whether  a  dynamic  uniaxial  stress-strain  relation  can  be  derived  from 
isothermal  or  adiabatic  quasi-static  conditions,  but  can  the  thermodynamic  condi¬ 
tions  represented  by  the  directly  obtained  relation  be  identified?  The  significance 
of  temperature  as  a  parameter  of  the  material  constitutive  relation  must,  then, 
be  assessed. 


It  has  been  recognized  that  the  elastic-plastic  stress  wave  represents 
an  irreversible  thermodynamic  process.  The  Hugoniot  relations  imply  that  the  work 
of  dynamically  stressing  a  material  from  aQ  to  is  given  by  the  area  under  the 
shock  chord  of  the  stress-strain  curve,  which  connects  the  two  stress  states 
(refer  to  Figure  D-16  on  page  D-26).  Note  that  this  quantity  differs  from  the 
area  under  the  stress-strain  curve  itself,  which  applies  for  gradual  loading.  Thus, 
the  plastic  shock  introduces  irreversible  dissipation  of  additional  mechanical  woik. 
It  has  been  shown  that  more  than  85  percent  of  the  work  goes  into  heat  dissipation, 
and  only  about  15  percent  into  permanent  changes  in  the  microstructure.  Con¬ 
sequently,  it  would  appear  that  the  final  str...,  state,  0^ ,  should  not  lie  on  a 
continuously  loaded  adia;  .tic  curve,  but  on  a  curve  modified  by  the  additional 
heat  release. 


(E-2) 

Leev  '  has  elaborated  further  on  the  above  by  proposing  that  the 
extremely  high  pressures  and  finite  dilatations  (up  to  25  percent  strain  at  400  kbar 
for  aluminum),  present  in  one-dimensional  strain  experiments,  introduce  dominant 
thermomechanical  coupling  effects  making  temperature  an  essential  variable  in  the 
equation-of-state.  If  this  assessment  applies  to  the  stress-strain  regitr  ;  of 
interest  in  the  study  of  spall,  the  accuracy  of  stress  predictions  from  rear 
surface  motion  measurements  would  be  lessened  considerably. 
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However,  other  investigators  have  given  consideration  to  the  thermo¬ 
dynamic  aspects  of  the  wave  propagation  and  have  developed  convincing  arguments 
in  favor  of  neglecting  temperature  effects  at  the  lower  stress  levels  of  the 
elastic-plastic  regime.  It  will  be  recalled  that  even  for  the  hydrodynamic  shock 
regime,  close  approximations  of  material  behavior  were  possible  by  neglecting 
temperature  effects  in  the  establishing  of  the  Hugoniot  material  (Appendix  C). 


Considerations,  similar  to  those  in  the  foregoing,  have  been  extended 
to  the  elastic-plastic  case.  The  basic  arguments,  developed  by  Horland,  are 
closely  paralleled  in  the  following  presentation.  Let  the  irreversible  energy 
dissipated  per  unit  mass  in  a  small  compression,  de,  be  dQ.  Conservation  of  work 
and  energy  gives 

dQ  -  dE  -  dW;  (E-2) 

where  dE  is  the  increase  in  internal  energy,  and  dW  is  the  work  done  by  the  applied 

/  T-t  f  V 

stress,  Taylor  and  Farren'1  '  have  shown  that  less  than  15  percent  of  the  plastic 
work  goes  into  producing  microscopic  structural  changes  in  the  material,  so  that 
heat  dissipation  per  unit  mass  represents  more  than  85  percent  of  dQ. 


On  integrating  over  the  full  compression,  the  increase  of  E  is  established 
from  Equation  D-61  as 

E1  "  Eo  “  7p^  <al  +  ao)  (ei  ’  eo)  •  02-3) 

The  work  increment  per  unit  mass  is 


dW 


(E-4) 


Hence,  the  total  applied  work,  ignoring  the  factor  1/p,  is  the  area  under  the 
corresponding  stress-strain  curve.  The  area  under  the  shock  chord  represents  the 
increase  in  E;  the  difference  represents  Q,  the  energy  dissipated. 


Referring  again  to  Figure  D-16,  it  is  evident  that  the  area  between  the 
stress-strain  curve  and  the  shock  chord  will,  in  general,  be  quite  small  if  the 
stress  jump  is  small  and  the  upward  curvature  of  the  hydrostat  is  negligible,  or 
if  work-hardening  effects  are  absent  or  linear  with  increasing  load.  An  estimate 
of  Q  may  be  found  by  assuming  that  the  stress-strain  relation  through  the  plastic 
shock  is  the  adiabatic  stress-strain  curve.  A  stress-strain  relation  of  the 
following  form  is  assumed: 

—  -a3  il+dje+dje  +  d3e  +  ...)  ; 


(E-5) 
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where  a,  f),  d^  ,  62  ...  are  determined  from  Bridgeman's  formula  fitting  hydro¬ 
static  material  data.  Integrating  the  above  between  the  yield  point,  o^,  and  the 
arbitrary  upper  stress  level  gives 

a  -  oy  -  3  a  3  [e  +  |  dxe2  +  j  d2e3  + . ]*  .  (K-6) 

Then , 


+  Oo)  -  {aY  -  3  «  3  fey  +  \  d^2  +  *“]}+  3  “  P  {|  +  eo)  +  \  dx  + 

eo2)  +  — }  •  (E-7) 


Also, 


J  lad«  -  (ex  -  eG)-[ay  -  3  a  3  [sY  +  \  d\  eY2  +  — -l}  +  3  “  P  {3  -  e02 )  +  3  di 

(el3  “  eo3)  +  —  j-  •  (E-8) 


Neglecting  terms  of  e  or  higher  gives 


(e  _e  ,3 

*  3p.  12  (  1  c°}  * 


Using  reasonable  values  for  aluminum,  Morland  establishes  that  for  Gj  -gq  <.  0.05, 
the  shock  change  is  virtually  adiabatic  (i.e.,  dQ  *  0),  so  that  the  material  behind 
the  shock  will  satisfy  the  original  adiabatic  stress-strain  relation. 


The  thermodynamic  considerations,  reviewed  in  the  foregoing,  establish 
the  validity  of  assuming  that  the  elastic-plastic  wave  process  is  adiabatic.  How¬ 
ever,  the  importance  of  temperature  on  the  wave  propagation  has  yet  to  be  assessed. 
Thus,  the  fact  that  the  process  is  essentially  adiabatic  indicates  that  a  tempera¬ 
ture  rise  occurs  under  the  action  of  the  shock  front.  The  resulting  interaction 
of  the  temperature  on  the  stress  field--the  couoled  thermomechanical  problem--has 
not  yet  been  solved  for  the  plastic  wave  case.  In  fact,  only  slight  progress  has 
been  made  in  the  case  of  thermoelastic  wave  propagation.  It  is  not  known,  for 
example,  whether  thermo-plastic  effects  might  be  confused  for  strain-rate  effects 
in  materials.  The  general  procedure  has  been  to  assume  that  the  propagating 
stress  field  is  unaffected  by  the  temperature  changes. 

Recently,  Lee  and  Liuv  '  and  Lee  and  Wierzbicki'  '  have  generalized 
elastic-plastic  theory  to  include  finite  strain  and  temperature  effects.  Although 
only  the  governing  equations  were  developed,  with  numerical  solutions  yet  to  be 
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obtained,  application  of  their  theory  to  specific  stress-wave  problems  may 
soon  yield  a  more  definite  assessment  of  the  importance  of  temperature  on  the 
elastic-plastic  wave  problem. 

Work-Hardening  and  Bauschinger  Effects 

As  was  evident  in  the  review  of  elastic-plastic  wave  propagation  in 
Appendix  D,  the  structure  of  the  propagating  stress  wave  is  a  function  of  the 
material  stress-strain  relation.  This  relation  will  be  characterized,  in  the 
absence  of  strain-rate  effects,  by  the  nature  of  its  loading  and  unloading.  Does 
the  material  have  a  sharp  or  ill-defined  yield  point?  Does  work  hardening  occur? 

Is  curvature  of  the  hydrostat  significant?  Is  a  Bauschinger  effect  present? 

If,  b'*  experimental  methods,  the  precise  shape  of  the  stress-strain  carve 
is  known,  no  assessment  of  the  above  effects  is  required  to  complete  the  analytical 
process  of  computing  the  interior  stress  field  in  a  plate-impact  t*st.  Yet,  in  the 
earlier  study  of  the  wave  propagation  and  interior  stress  field  'Appendix  D, 
pages  D-33-D-45) .even  the  use  of  a  simple  material  model  (elastic,  perfectly-plastic) 
led  to  a  complicated  wave  structure  due  to  many  wave  interactions.  However,  the 
availability  of  computer  codes  for  analyzing  the  stress-wave  propagation,  in¬ 
cluding  all  of  the  wave  front  and  interface  interactions,  removes  the  analytical 
difficulties  resulting  from  complexities  in  the  material  stress-strain  relation. 
Hence,  the  presence  or  absence  of  the  aforementioned  effects,  which  serve  to 
distinguish  various  materials,  presents  no  fundamental  difficulty  in  solving  the 
stress-wave  problem. 

The  preceding  remarks  are  presupposed  on  the  basis  of  knowing  the  material 
stress-strain  relation.  However,  if  inexactitude  exists  in  the  experimental  methods 
for  establishing  the  stress-strain  relations,  there  will  be  a  consequent  uncer¬ 
tainty  as  to  the  interior  stress  field.  The  present  method  for  establishing  the 
stress-strain  relation  (reviewed  on  pages  D-33-D-45)  is  to  measure  the  rear  sur¬ 
face  motion  in  a  plate-slap  experiment  and,  by  trial-and-error  solution  of  the 
stress-wave  problem,  to  deduce  a  stress-strain  curve  that  will  yield  a  theoretical 
prediction  of  the  rear  surface  motion  in  accord  with  the  experiment.  It,  there¬ 
fore,  becomes  important  to  assess  the  accuracy  of  the  stress-strain  relation 
obtained  from  the  plate-slap  experiment  itself. 
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The  first  remark  to  be  made  in  this  context  concerns  the  sensitivity  of 

the  stress-strain  relation  to  erroneous  measurement  of  surface  motion.  In  this 

(E-7) 

regard,  it  has  been  pointed  out  by  Butcher  and  Canon'  '  that  errors  in  the 
measurement  of  surface  velocity,  when  converted  to  material  velocity,  cause  data 
points  to  move  along  the  dynamic  stress-strain  curve  rather  than,  say,  parallel 
to  the  stress  or  strain  axis.  This  means  that  errors  in  surface  motion  measurement 
are  reflected  as  much  smaller  errors  in  the  determination  of  the  stress-strain 
curve. 

Further  assessment  may  be  divided  into  consideration  of  the  experimentally 
obtained  loading  and  unloading  information  about  the  material.  There  is,  in  fact, 
basis  for  confidence  in  the  rear  surface  motion  measurement  approach  to  establishing 
the  loading  characteristics  of  the  stress-strain  curve.  Thus,  the  magnitude 
of  the  elastic  precursor  establishes  the  material  yield  point.  The  propagation 
velocity  of  the  plastic  front  establishes  the  shock  chord  on  the  stress-strain 
curve,  and,  by  performing  several  exper_ments  at  differing  plate-impact  velocities, 
several  stress-strain  points  on  the  loading  curve  may  be  determined.  The  elastic- 
plastic  compression  waves  that  establish  the  above  data  are  the  first  waves  to 
arrive  at  the  free  surface,  and  they  have  done  so  with  a  minimum  of  self- 
interaction  (see  pages  D-33-D-45  for  the  discussion  of  the  interactions  that 
have  occurred) . 

However,  the  deduction  of  the  unloading  characteristics  of  a  material 

from  rear  surface  measurements  is  less  certain.  Thus,  the  unloading  information 

appearing  in  the  form  of  rear  surface  motion  is  due  to  waves  reflecting  from  the 

back  surface  of  the  projectile-target  combination.  These  waves  arrive  at  the 

rear  surface  only  after  traveling  as  compression,  then  release,  and  finally  as 

(E-l) 

tension,  undergoing  numerous  interactions  in  the  process.  Barker,  et  al, 
discuss  the  necessity  for  performing  tests  using  various  target-projectile 
thickness  ratios.  Use  of  various  ratios  would  help  to  establish  unique  unloading 
information,  which  is  difficult  to  acquire  because  of  multiple  reflections  and 
interactions.  The  necessity  for  performing  multiple  experiments,  and  resultant 
indirectness  in  establishing  the  scress-strain  curve,  casts  doubt  on  the  accuracy 
of  unloading  and  reverse  loading  information  established  by  rear  surface  motion 


measurements . 
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Until  recently,  the  primary  method  for  determining  material  properties 
was  by  inference  from  rear  surface  motion  measurements.  Refinements  in  measuring 
technique,  utilizing  laser  velocity  interferometers,  have  been  reported 
[Karnes''  ;  ] .  However,  Barker'  has  recently  reported  on  measurements  of 
material  particle  velocities  in  the  interior  of  a  shocked  solid  wherein  the  laser 
interferometer  was  used  in  conjunction  with  a  transparent  quartz  window  on  the 
specimen  to  obtain  the  interior  data.  < ransparent  window,  of  similar  mechanical 

impedance  as  the  shocked  material  (ale  s  attached  to  the  specimen.  The 

interface  motion  was  obtained  with  the  ini..  meter  and  was  assumed  to  depict  the 
motion  of  a  truly  interior  point  of  the  material.  Although  the  complete  validity 
of  the  technique  has  not  been  established,  since  questions  of  slight  impedance 
mismatch  effects  and  the  light  transmission  characteristics  of  quartz  under  dynamic 
loads  remain,  it  does,  apparently,  enable  the  elastic-plastic  wave  structures  to 
be  established  directly,  before  wave  reflections  and  interactions  have  occurred. 

As  a  consequence  of  this  development,  direct  information  on  the  unloading  charac¬ 
teristics  of  materials  is  now  available. 

In  Figure  E-l,  the  data  obtained  by  Barker  on  the  particle  velocity¬ 
time  history  at  an  interior  surface  is  shown  for  three  different  tests.  Consider, 
in  the  figure,  the  general  nature  of  the  loading  portion  of  the  wave  and  of  the 
unloading  portion.  It  is  to  be  noted  that  at  the  higher  stress  level  shots 
(Nos.  926  and  927),  the  plastic  front  is  discontinuous,  but  the  low  stress  level 
of  Shot  922  yields  a  rather  smooth  stress  variation.  This  is  cited  by  Barker  as 
evidence  cf  a  strain-rate  effect  and  is  discussed  later.  The  unloading  portion 
of  the  wave,  not  heretofore  directly  measured,  exhibits  considerable  lack  of 
structure.  This  also  has  been  partially  attributed  to  a  strain-rate  effect.  How¬ 
ever,  comparing  predicted  and  observed  results,  it  is  evident  that  considerably 
more  release  wave  structure  was  predicted  by  some  material  models. 

In  order  to  obtain  comparable  experimental  and  theoretical  resuits, 

Barker  found  it  necessary  to  postulate  a  pronounced  Bauschinger  effect  for  aluminum. 
Two  other  material  models  did  not  include  this  effect.  The  simplified  model, 
identified  as  Y  *  2.4  kbar  model  in  Figure  E-l  (i.e,,  a  constant  yield  strength  of 
2.4  kbar),  postulated  no  Bauschinger  effect  and  no  work  hardening.  Although  the 
loading  behavior  predicted  by  this  model  is  as  satisfactory  as  the  more  complicated 
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models,  considerable  deviation  from  the  observed  data  is  predicted  on  unloading. 

A  second  model,  without  a  Bauschinger  effect  but  including  work  hardening 

(Y  «  2.4  +  0.05P,  P  =  hydrostatic  stress),  also  has  significant  discrepancies  in 

the  predicted  unloading. 

In  order  that  the  o^^icted  release  wave  take  on  the  observed  form,  it 
was  necessary  for  the  rev  -  yield  strength  on  unloading  to  be  initially  small, 
then  to  increase-~at  first  rapidly  and  later  more  slowly.  This  amounted  to  the 
pronounced  Bauschinger  effect.  The  mathematical  expression  for  this  effect  was 
given  by  describing  the  ten  ile  yield  strength,  Y-p,  as  a  function  of  plastic 
strain,  e p,  according  to  the  form 

YT  «  4[1  -  exp  (-800  ep)].  (E-9) 

In  addition,  a  work-hardening  model,  i lightly  different  than  the  one  previously 
given,  was  used  for  the  loading  behavior.  It  was  assumed  that  the  compressive 
yield  strength  increased  from  2.4  kbar  to  3.0  kbar.  The  resulting  wave  predic¬ 
tion  is  shown  as  the  dark  iine  in  Figure  E-l.  It  is  seen  that  considerably  better 
agreement  between  truury  and  test  exists.  Of  "course,  the  fine  structure  of  the 
predicted  wave  shape  is  a  result  of  the  numerical  analysis  which  uses  the  method 
of  characteristics.  This  breaks  a  stress  pulse  into  stress  increments  or  steps, 
as  discussed  in  earlier  sections  of  this  report,  and  of  necessity  leads  to  a  pre¬ 
dicted  wave  shape  possessing  definite  structure.  The  three  stress-scrain  models 
used  are  compared  in  Figure  E-2. 

By  observing  Figures  E-2  and  E-l,  it  is  possible  to  assess,  at  least 
qualitatively,  the  effect  of  uncertainty  in  the  stress-strain  curve  on  stress 
pulse  prediction.  It  is  seen  that  all  cf  the  models  accurately  predict  the  maximum 
stresses,  aid  that  at  the  highest  stress  level  (Shot  926  of  Figure  E-l),  all  of 
the  models  adequately  predict  the  loading  behavior  in  general.  At  the  lower 
stresses  (particularly  in  Shot  922),  discrepancies  in  the  plastic  shock  front 
occur,  and  these  are  thought  to  hd  due  to  strain-rate  effects.  Thus,  at  lowest 
stress  level,  tone  of  the  mrdels  properly  depict  the  nature  of  th-  plastic  com¬ 
pression  front. 

Considering  the  unloading  portions  of  the  various  waves,  it  i?  jteen  that, 
significant  discrepancies  occur  in  the  material  models  that  do  not  contain 
Bauschinger  effects.  Thus,  fairly  significant  differences  occur  between  measured 


FIGURE  E-2.  LOADING  AND  UNLOADING  STRESS-STRAIN  PATHS 
FOR  SHOT  927  ACCORDING  TO  THREE  MODELS  OF 
MATERIAL  BEHAVIOR^2-9) 
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and  predicted  unloading  at  all  three  stress  levels  for  the  constant  yield  model 
(Y  *  2.4  kbar)  and  t!  work-hardening  model  (Y  >*  2.4  +  0.05  P) .  The  inclusion 
of  the  Bauschinger  effect,  in  the  manner  indicated  by  Equation  E-9,  yields  close 
agreement  between  the  predicted  and  measured  results.  Yet,  referring  to  Figure 
E-2,  it  is  apparent  that  the  percentage  differences  among  all  three  models  is 
quite  small.  This  suggests  that  the  stress  history  within  the  material, 
particularly  the  unloading  history,  is  a  sensitive  function  of  the  stress-strain 
relation.  Sint.  he  resultant  spall-producing  tensile  wavs  is  a  result  of  inter¬ 
action  between  unloading  waves,  it  becomes  important  to  have  an  accurate  knowledge 
of  the  unloading  behavior  of  a  material  to  be  able  to  predict  spall. 
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£  ‘.raln-Rate  Effects 

In  the  p-^vious  sections,  some  statements,  at  least  qualitative,  were 
possible  on  Rauschinger  or  work-hardening  effects  and  materiel  response.  In  fact, 
rather  than  actually  assessing  the  effects,  it  was  a  matter  of  establishing  the 
presence  or  absence  of  the  effects  by  experimental  methods  and  incorporating  them 
into  the  resultant  material  constitutive  relation.  The  situation  is  much  less 
certain  in  the  case  of  strain-rate  effects.  At  present  there  is  considerable 
controversy  regarding  the  degree  of  strain-rate  sensitivity  of  materials.  More 
basic  than  this  is  the  controversy  aroused  by  the  question  of  whether  strain 
rate  is  even  a  parameter  of  the  constitutive  relation  for  some  materials. 

Nevertheless,  it  is  a  fact  that  wave- propagation  experiments  in  rods  and 

thick  plates  have  yielded  results,  all  details  of  which  cannot  be  explained  by  a 

strain-rate  independent  theory.  For  example,  in  the  review  of  elastic-plastic 

wave  propagation  (Appendix  D) ,  particularly  the  part  referring  to  Figure  D-25, 

it  is  indicated  that  sharply  discontinuous  plastic  shock  fronts  are  expected  from 

(E-l) 

a  strain-rate -independent  theory.  However,  Barker,  et  al,  have  indicated 

that  their  experimental  results  were  a  consequence  of  stress  waves  having  the 
general  configuration  shown  in  Tigure  E-3.  Worthy  of  note  in  Figure  E-3  is  the 
lack  of  a  sharp  plastic  front,  with  the  strain  gradually  attaining  its  maximum 
value  as  though  a  relaxation  process  were  operative.  Referring  also  to  Figure  E-l, 
which  showed  Barker's  results,  it  is  evident  that  sharp  discontinuities  are 
absent,  particularly  at  the  lower  loadings. 


FIGURE  E-3.  QUALITATIVE  DRAWING  OF  A  WAVE  SHAPE  WHICH  INDUCED 
THE  OBSERVED  FREE  SURFACE  MOTION^-1' 

rhe  presence  of  a  strain-rate  effect,  in  terms  of  a  general  viscoplas tic 
relaxation  mechanism  would  tend  to  explain  qualitatively  the  results  obtained. 

This  may  be  illustrated  \  considering  a  material  that  exhibits  the  characteristics 
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of  a  simple  Voigt  viscoeiastic  model,  shown  in  Figure  E-4.  Thus,  it  is  seen  that 
the  application  of  a  sudden,  step- function-type  stress  to  the  material  would  re¬ 
sult  in  a  delayed  strain  response,  with  the  strain  approaching  its  final  value 
only  gradually.  CJ 


■fc- 

(0) 


(b) 


OF  237 


FIGURE  E-4  (a) .  SIMPLE  VOIGl’  VISCOELASTIC  MODEL 
(b) .  STRAIN  RESPONSE  TO  STEP  LOADING 

Although  certain  experimental  results  indicate  strain-rate  sensitivity 
in  materials,  and  the  operation  of  certain  viscoplastic  mechanisms  could  account 
for  some  of  the  observations ,  considerable  controversy  still  attaches  to  this 
subject.  In  fact,  interest  in  this  matter  arose  with  the  earliest  developments  of 
elastic-plastic  th-^ry  by  von  Karman,  ^E-1^  Taylor,  ^E~^  and  Rakhmatulin . 'E 
Thus,  these  earlie-  *„  developments  pertained  to  a  strain-rate-independent  theory 
of  elas tic-plastic  waves  ir.  rods  and  succeeded  in  predicting  the  essential  features 
of  such  wave  propagation,  e.g.,  the  existence  and  magnitude  of  an  elastic  precursor 
and  the  existence  of  a  plateau  of  constant  strain  behind  the  wave  front.  The 
theory  also  predicted  that  superimposed  stress  increments  would  propagate  at  the 
plastic-wave  velocity  given  by  the  local  tangent  modulus  of  the  stress-strain  curve. 
However,  experiments  by  Bell^  ^  and  others  showed  the  increments  to  propagate  at 
the  elastic  velocity.  In  addition,  the  theory  failed  to  account  f  the  increase 
in  yield  strength  noted  for  some  materials,  such  as  iron,  under  dynamic  loading. 

In  an  effort  to  remove  the  apparent  anomalies  of  the  theory, 

Malvern ^ ’  ^E  ^  developed  a  strain-rate-dependent  theory  where  the  basic 
constitutive  equation  took  the  form, 

Ee  *  a  +  g(o; f.)  ;  (E-10) 

wl  're  e  is  the  total  strain,  a  is  stress,  and  E  is  Young’s  modulus.  The  function 
g(a,e)  is  related  to  the  plastic  strain  rate  6p  by 

Etp  «  g(a,e)  .  (E-ll) 

Although  some  of  the  anomalies  w^re  removed  by  this  theory,  the  constant  strain 
plateau  predicted  by  the  rate-independent  theory  and  observed  experimentally  was  not 
predicted  by  the  rate-dependent  theory. 
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In  tne  ensuing  years,  there  have  apparently  developed  "schools"  of 
thinking  regarding  the  sensitivity  of  materials  to  strain  rate.  Bell  and  his 
co-workers  have  been  the  principal  proponents  of  the  strain-rate-independent 
behavior  of  materials;  Lindholm,^  ^  Malvern,  and  others  have  suggested  that 
materials  are  strain-rate  sensitive  to  one  degree  or  another.  Unfortunately, 
the  experimental  techniques  used  by  the  various  investigators  to  arrive  at  their 
differing  conclusions  on  dynamic  materials  properties  have  been  different.  Thus, 
the  difficulty  of  assessing  the  importance  of  strain-rate  effects  is  increased. 
Bell,  for  example,  has  utilized  his  unique  diffraction  grating  technique  for 
studying  elastic-plastic  wave  propagation  in  long  rods.  Lir.dholm  and  others  of 
the  rate-dependent  school  have  used  the  split  Hopkinson  bar  method.  Although 
members  of  each  school  have  developed  impressive  arguments  tending  to  show  the 
experimental  methods  of  the  other  school  at  fault,  a  definitive  set  of  compara¬ 
tive  experiments  has  not  been  conducted  to  settle  this  mattar. 

It  has  been  suggested  that,  in  fact,  both  schools  may  be  correct--under 
the  proper  conditions.  Dorn  and  Hauser^  17^  have  postulated  that,  under  proper 
conditions,  certain  materials  may  be  essentially  strain-rate  independent;  and 
other  materials,  under  proper  conditions,  may  be  quite  sensitive  to  rate  effects. 
Of  course,  there  exist  intermediate  regions  of  varying  degrees  of  sensitivity. 

The  factors  differentiating  the  degrees  of  sensitivity  were  postulated  to  be  the 
dislocation  activation  mechanisms.  Thus,  strain-rate-sensitive  materials  were 
postulated  to  be  those  whose  dislocation  mechanisms  of  deformation  could  be 
thermally  activated,  whereas  rate  insensitive  materials  had  mechanisms  that  were 
athermally  or  stress  activated. 

Although  the  above  postulation  of  intermediate  position  between  the 
strain-rate  schools  certainly  appears  reasonable,  the  problem  of  establishing  the 
physical  classifications  of  materials,  including  the  various  subtle  transition 
regions,  has  not  yet  been  solved. 

Accepting  the  existing  uncertain  state  of  affairs  regarding  strain-rate 
sensitivity  of  materials,  it  is  instructive  to  consider  the  recent  work  of  Butcher 
and  Kames^"^  on  strain-rate  effects  in  metals.  In  their  work,  the  strain-rate 
theory  of  Malvern  was  applied  to  the  one-dimensional  strain  configuration  of  the 
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plate-slap  experiment.  Predicted  stress-strain  curves  were  determined  and  found  to 
be  more  in  accord  with  experimental  results  than  rate- independent  theory.  Thus, 
the  Malvern  theory,  which  yielded  some  anomalous  results  in  the  one-dimensional 
stress  application  (no  constant  strain  plateau),  successfully  predicted  the  es¬ 
sential  features  of  the  experimentally  (bserved  response  in  the  one-dimensional 
strain  situation.  However,  it  was  also  pointed  out  by  Butcher  and  Karnes  that 
certain  experimental  factors,  in  particular  a  slight  lack  of  planarity  between  the 
impacting  plates,  were  capable  of  producing  results  resembling  strain-rate  effects. 
It  was  also  indicated  the-  increased  distance  of  travel  tended  to  decrease  rate 
effects  in  the  wave  front,  permitting  the  wave  to  be  treated  as  essentially 
strain-rate  independent.  Although  future  study  of  strain-rate  effects  utilizing 
the  plate-slap  configuration  will  certainly  be  done,  the  general  conclusions 
reached  by  Butcher  and  Karnes  appeared  to  be  that  strain-rate  effects  were  rather 
slight  under  plate-slap  conditions  and  that,  furthermore,  the  presence  or  plate 
tilt  and  propagation  distance  effects  was  capable  of  obscuring  rate  effects. 

The  basis  for  the  relative  insensitivity  of  the  one-dimensional  strain 

configuration  (in  comparison  to  the  one-dimensional  stress  situation)  to  strain- 

(E-l) 

rate  effects  has  been  discussed  by  Barker,  et  al.  In  this  configuration,  the 

majority  of  the  strain  is  diletational  and  not  deviatoric.  However,  it  is  the 
deviatoric,  or  shearing-type  deformation,  .hat  causes  dislocation  motion  and  the 
consequent  strain-rate  effect.  It  follows  that  strain-rate  effects  will  be  con¬ 
siderably  less  significant  in  such  configurations, 

Despit  *  the  considerable  research  that  has  been  directed  toward  as¬ 
certaining  the  sensitivity  of  materials  to  rate  effects,  it  appears  that  material 
constitutive  relations,  including  these  effects  and  accounting  for  the  various 
experimental  observations,  have  not  been  established  although  many  have  been 
postulated.  Thus,  in  contrast  to  the  previously  discussed  Bauschinger  and  work¬ 
hardening  effects,  where  explicit  inclusion  of  these  effects  in  the  stress-strain 
relation  became  possible  on  the  basis  of  experimental  observations,  similar 
developments  have  not  been  made  in  regard  to  rate  effects. 

Nevertheless,  when  one-dimensional  strain  considerations  govern,  it 
has  been  indicated  in  the  previous  discussion  that  strain-rate  effects  appear 
to  be  of  little  significance  in  predicting  the  essential  features  of  the  wave 
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propagation.  Although,  as  is  evident  from  Barker's  work,  some  detail  of  the 
wave  structure  is  lost,  the  essential  features  of  the  waves,  such  as  stress 
amplitude,  are  predictable  by  a  strain-rate-independent  theory. 

Finite  Strain  Effects 


The  final  topic  of  discussion  regarding  various  subsidiary  effects 
pertains  to  the  kinematic  description  of  the  deformation.  The  theoretical  develop¬ 
ment  of  elastic-plastic  theory  reviewed  in  this  report  utilized  a  description  of 

du 

the  Lagrangian  strain  as  ex  =  /dx.  Such  a  representation  is  tlu.  infinitesimal 
strain  definition  and  neglects  quadratic  terms  in  the  displacement  gradients.  Yet, 
under  plastic  deformation,  large  strains  may  exist  and,  thus,  may  require  the  more 
accurate  kinematical  description  of  strain.  Lee,  '  for  example,  has  suggested 
the  necessity  of  including  finite  strain  effects  in  the  theory,  since  strains  of 
25  percent  are  evidently  obtained  in  some  tests.  The  basis  for  using  the  infini¬ 
tesimal  definition  should  be  reviewed. 

If,  in  fact,  strains  of  the  magnitude  indicated  above  were  typical  of 
those  attained  in  che  low  pressure  regions  of  interest  in  elastic-plastic  theory, 
the  question  of  finite  strain  effects  would  need  to  be  examined  with  considerable 
care.  However,  such  large  strains  are  associated  more  often  with  pure  hydro- 
dynamic  shock-type  loading  where  extremely  high  pressures  are  attained,  whereas, 
the  stress  regions  of  interest  in  elastic-plastic  wave  propagation  are  of  the 
order  10-50  kbar,  and  the  associated  strains  are  2-4  percent.  Under  these  con¬ 
ditions,  the  squares  of  the  displacement  gradient  are  negligible  compared  to  the 
displacement  gradient.  It  might  be  noted  further  that,  under  strains  greatly  in 
excess  of  the  2-4  percent  level,  the  differences  ir.  the  principal  stresses  becomes 
quite  small  compared  to  the  stresses  themselves.  Consequently,  the  elastic-plastic 
theory  becomes  little  more  than  a  small,  superimposed  correction  to  a  hydrodynamic 
model  of  the  material.  Hence,  where  it  becomes  necessary  to  include  the  more 
complicated  finite  strain  effects,  a  corresponding  simplification  may  be  made  in 
the  nature  of  the  constitutive  relation. 

On  the  basis  of  qualitative  arguments,  it  appears  safe  to  neglect  finite 
strain  effects  on  the  elastic-pla.  «■*  c  wave  propagation.  It  is  instructive  to 
exp’ore  this  topic  further  to  see  if  a  more  precise  assessment  may  be  attached 
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to  finite  strain  effects.  It  3nould  first  be  noted  that,  in  general,  when  the 
deformation  of  a  body  becomes  sufficiently  large,  it  can  no  longer  be  suitably 
characterized  in  terms  of  either  "engineering"  or  "true"  strain*  This  is  due  to 
the  fact  that  these  measures  of  deformation  are  not  inherently  tensor  quantities 
and  are  inappropriate  for  the  establishment  of  a  constitutive  relationship  that 
is  to  be  independent  of  the  choice  of  coordinates.  In  addition,  if  Lagrangian 
coordinates  are  being  used,  the  Eulerian  stress  tensor  of  linear  elasticity  theory 
must  be  distinguished  from  and  replaced  by  a  Lagrangian  or  Kirchhoff  stress  tensor 
based  on  predeformation  geometry. 

For  the  case  of  wave  propagation  in  uniaxial  strain,  there  is  no  need  to 
distinguish  between  Eulerian  and  Lagrangian  normal-to-surface  stresses.  Also,  as 
far  as  stress  and  displacement  wave  propagation  is  concerned,  it  makes  no  difference 
whether  the  dynamic  stress-strain  relationship  is  given  in  terms  of  the  components 
of  a  finite  strain  tensor  or  engineering  strain. 


To  see  this,  consider  a  material  with  dynamic  stress  versus  finite-strain 
curve  defined  by 


<?x  “  ax(E)  ; 
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where 


ax  ■  stress  normal  to  the  direction  of  wave  propagation, 

E  »  finite  strain  normal  to  the  direction  of  wave  propagation  (the 
only  nonzero  component  of  the  strain  tensor). 

Suppose,  also,  that  the  finite-strain  component,  E,  and  the  engineering  strain 
component,  ex,  are  related  by 

E  »  E(ex)  ; 

where  ex  is  the  previously  defined  infinitesimal  strain,  given  by 


(E- 1 3  > 
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u  *  particle  displacement  in  the  direction  of  wave  propagation, 
x  *  Langrangian  coordinate  in  the  direction  of  wave  propagation. 
The  expressions  for  conservation  of  mass  and  momentum  are,  respectively, 


Po 


p(l  +  ex) , 


(E-15) 
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Nc~f  assuming  continuous  first  derivatives, 


3x  Be  Bx  Be  3ex  Bx 


(E- 16) 


However , 


Be 

Be. 


30^ 

Bex 


so  that 

2 

Bax  _  Box  Bex  Bcx  B  u 
Bx  3ex  Bx  3ex  ^  2 
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Consequently,  the  governing  equation  reduces  to  that  which  would  be  obtained  if  only 
engineering  strain  were  considered,  i.e., 


3t 


2 
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where 
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c  = 


Po  Se 
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Thus,  while  the  stress  versus  finite-strain  and  stress  versus  engineering-strain 
relationships  may  be  of  entirely  different  form,  the  spatial  derivative  of  ax  and, 
hence,  the  governing  equation  of  motion  is  unchanged. 


In  order  to  illustrate  the  above  comments  further,  it  will  be  instructive 
to  consider  wave  propagation  in  two  rather  ideal  materials,  ^or  the  purpose  of 
this  discussion,  the  Lagrangian  strain  tensor  will  be  used.  The  only  nonvanishing 
component  of  this  tensor  will  be  E,  the  strain  normal  to  the  direction  of  wave 
propagation.  For  uniaxial  strain, 


E 

First,  consider  a  material 
strain,  ex;  i.e., 


\  2 


=  e 


,1  2  3u  1  /  3u\ 

+  o  ex  =  •ET  +  7  [aZ J  • 
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hat  is  linearly  elastic  with  respect  to  engineering 
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ax  =  K*ex  ,  (E -21) 

where  K*  is  the  effective  dynamic  modulus  in  uniaxial  strain.  Expressing  ex  in 
terms  of  E  gives 


ex  =  +  J 1  +2E  -  1,  (ex  >  -1,  E  >  -1/2). 
The  stress  versus  finite-strain  relationship  becomes 


(E- 22) 


°x  =  K*  {+  V1  +  2E  -  0- 


(E-23) 
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Now 


and 


Hence, 
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Thus,  the  equation  of  motion  is  the  same  as  it  would  have  been  * f  finite  strain 
were  not  considered,  i.e., 


dji 

at2 


2  a  u  2  K* 


ax 


2' 


Po 
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Now  consider  a  material  that  is  linearly  elastic  with  respect  to  finite 
strain.  Rubber-likc  materials  can  often  be  described  in  this  way.  Let  the  stress 
versus  finite-strain  relationship  be  given  by 

0x  m  K  E  ,  (E-28) 

where  K  is  an  effective  dynamic  modulus  in  uniaxial  strain.  Now 


*2*  .  |F  ^ 

ax  *  3x 


and  the  equation  of  motion  becomes 


K  (1  +  ex) 


(E-29) 


where 


c2  -  K,lLt.eX?..  .  K 

Po  P  ' 


(E-30) 
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The  interesting  and  perhaps  surprising  result  is  that,  in  this  type  of 
material,  large  tensile  disturbances  propagate  at  a  higher  velocity  than  small 
tensile  disturbances,  and  large  compressive  disturbances  propagate  at  a  lower 
velocity  than  small  compressive  disturbances.  Such  behavior  has  even  been 


E-20 


observed  experimentally.  Kolsky, ( E— 19)  for  example,  reports  that  the  propagation 
velocity  of  small  tensile  pulses  in  neoprene  filaments  increased  by  a  factor  of  ten 
who”  the  filament  was  prestretched  to  six  times  its  original  length. 


This  behavior  is  not  so  surprising  if  it  is  realized  that  a  material, 
linearly  elastic  with  respect  to  finite  strain,  is  not  linear  with  respect  to 
engineering  strain.  Using  the  Lagrangian  definition  of  finite  strain,  the  ctress- 
engineering-strain  curve  would  have  the  parabolic  form  shown  in  Figure  E-5.  The 
"sound"  speed  for  a  material  with  this  type  of  stress-engineering-strain  relation¬ 
ship  is  given  by 


c2  _  i_  _  K  (1  +  ex) 
Po  ^ex  Po 


(E-.',2) 


This  is  identical  to  the  sound  speed  obtained  previously  with  the  stress-finite 
strain  relationship  (Equation  E-31) . 


cr  =  KE  =  K(e  +  £e2) 
<r  =  Ke 
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FIGURE  E-5.  NONLINEAR  AND  ENGINEERING 
STRESS-STRAIN  CURVES 

From  the  discussion  above,  it  may  then  be  concluded  that  if  the  stress 
versus  engineering-strain  relationship  for  uniaxial  strain  is  known,  there  i3  no 
need  to  add  a  finite  strain  correction  for  large  values  of  strain  since  the 
equation  governing  displacements  would  be  unchanged.  The  magnitude  or  amplitude 
of  the  strain  pulse,  however,  does  change  according  to  how  strain  is  defined;  i.e., 
the  Lagrangian  strains  are  always  larger  than  the  engineering  strains.  But  these 
quantities  are  not  generally  of  interest.  It  should  perhaps  also  be  emphasized 
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that  the  magnitude  of  the  stress  pt_lse  does  not  change.  Since  the  use  of  a  stress 
finite-strain  relationship  was  shown  to  have  no  effect  on  the  equation  governing 
displacement  wave  propagation,  application  of  the  Rankine-llugoniot  jump  conditions 
to  determine  stress  from  shock  and  particle  velocity  measurements  would  be  no  less 
valid  than  it  was  when  only  engineering  strain  was  considered.  Thus,  although 
E  (x,t)  ^  ex  (x,t),  the  stress  calculated  as  a  function  of  E  and  the  stress 
calculated  as  a  function  of  ex  are  equal  because  of  the  different  moduli  involved. 
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